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INTRODUCTION 

This  research  project  involves  the  analysis  of  an  existing  clinical  set  (~150  patients)  of  frequency- 
domain  optical  mammograms  to  assess  the  performance  of  this  approach  to  breast  cancer  detection. 
The  analysis  of  the  breast  images  is  complemented  by  theoretical  and  experimental  studies  to 
characterize  the  proposed  algorithms  of  image  processing.  The  objective  of  this  research  is  to 
identify  the  strengths  and  the  weaknesses  of  the  current  instrument  design,  and  to  guide  the  design 
of  new  optical  instrumentation  for  breast  cancer  detection. 


BODY 

Approved  statement  of  work 

The  approved  statement  of  work  for  this  project  is  the  following: 

Task  1. 

Compute  and  analyze  the  edge-corrected  optical  mammograms  for  all  150  patients  (Months  1-18) 

a.  Install  the  optical  mammography  software  for  the  SUN  workstation; 

b.  Build  the  edge-effect-corrected  images  (iV-images)  from  the  amplitude  and  phase  images; 

c.  Carry  out  a  comparison  between  the  TV-images  and  the  tumor  diagnosis  for  each  patient; 

d.  Analyze  the  wavelength  dependence  of  the  N  parameter  and  the  ac  amplitude  for  different 
kinds  of  tumors; 

e.  Collect  a  printed  summary  of  all  the  2,400  N- images  (150  patients,  2  breasts,  2 
projections,  4  wavelengths)  and  the  spectral  features  of  all  the  tumors. 


Task  2. 

Perform  the  optical  measurements  on  breast-like  phantoms  (Months  19-24) 

a.  Prepare  the  breast-like  phantoms  (optical  inhomogeneities  +  strongly  scattering 
background); 

b.  Collect  and  analyze  the  optical  data; 

c.  Repeat  the  experiment  for  a  variety  of  sizes,  shapes,  optical  contrasts,  and  inhomogeneous 
backgrounds. 


Task  3. 

Apply  the  method  for  the  quantification  of  the  tumor  optical  properties  and  construct  the  spectral 
chart  for  benign  and  malignant  tumors  (Months  25-36) 

a.  Find  the  optical  properties,  and  the  hemoglobin-related  parameters  of  the  tumors; 

b.  Build  charts  of  benign  and  malignant  regions  in  the  iviV  plane  at  four  wavelengths; 

c.  Determine  whether  malignant  tumors  show  a  specific  spectral  signature. 

The  first  Grant  period  (year  01)  was  devoted  to  Task  1,  and  the  results  obtained  were  reported  in  the 
first  year  annual  report.  The  second  Grant  period  (year  02)  was  devoted  to  an  extension  of  Tasks  lb 
and  Id  and  to  Task  2. 
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Task  l.b. 

During  year  01,  we  have  obtained  the  so-called  N-images  for  the  whole  dataset  of  optical 
mammograms.  We  briefly  recall  that  the  N-images  are  the  result  of  an  algorithm  of  image 
processing  previously  developed  by  us  to  minimize  the  effect  of  the  breast  thickness  variability  on 
the  optical  mammograms.  The  resulting  images  display  a  dimensionless  parameter  (called  N)  that 
is  indicative  of  the  optical  density  of  the  breast.  We  have  found  that  N-images  enhance  the  tumor 
detectability  with  respect  to  amplitude  or  phase  images  (i.e.  with  respect  to  images  obtained  from 
the  raw  optical  data).  After  having  examined  the  results  for  the  whole  dataset,  we  have  found  that 
all  of  the  outstanding  image  regions  associated  with  cancer,  benign  tumors,  and  blood  vessels  are 
characterized  by  an  increased  value  of  N.  Such  an  increased  value  (peak)  is  always  associated  with 
a  negative  second  spatial  derivative.  On  the  contrary,  we  have  observed  that  areas  that  introduce 
edge  effects  (i.e.  the  artifacts  associated  with  the  breast  thickness  reduction  by  the  edge  of  the  breast 
image)  are  usually  associated  with  a  positive  second  spatial  derivative  of  N. 

As  a  result  of  the  above  findings,  during  year  02,  we  have  started  developing  a  new  scheme 
of  image  processing  based  on  the  second  derivative  of  N.  This  image  processing  approach  consists 
of  the  following  steps: 

I)  Smooth  the  original  N  image  by  a  low-pass  spatial  filter; 

II)  Calculate  the  average  of  the  second  spatial  derivative  along  eight  directions  through 
each  image  pixel  (N”); 

III)  Set  all  of  the  image  pixels  corresponding  to  a  positive  value  of  N”  to  white; 

IV)  Set  all  of  the  pixels  out  of  the  breast  to  black; 

V)  Display  the  image  pixels  with  N”<0  using  a  gray  scale  palette  where  white  corresponds 
to  N”=0  and  black  corresponds  to  the  maximum  value  of  N”. 

This  algorithm  of  data  processing  keeps  all  of  the  strengths  of  the  N-imaging  and  further  improves 
on  the  quality  of  the  image  display.  In  fact,  (a)  the  correction  of  the  edge  artifacts  afforded  by  the 
N-image  is  further  improved  by  taking  the  second  derivative  (because  the  edge  zone  will  be 
typically  set  to  white);  (b)  the  automatic  generation  of  N-images,  which  is  a  very  important  feature 
for  the  practical  examinations  of  the  images  in  a  clinical  setting,  is  also  featured  by  the  second 
derivative  images;  (c)  secondary  peaks  in  the  N-images,  which  may  be  buried  within  the  main  peak, 
appear  with  sharper  contrast  in  N”  images,  and  (d)  the  shape  of  the  areas  of  increased  N  are  much 
better  identified  in  the  N”  images  rather  than  the  N  images.  Figure  1  reports  representative  N 
images  of  the  optical  mammograms  at  690  nm  collected  on  a  patient  affected  by  cancer.  The 
cancer,  3  cm  in  size,  appears  as  a  blurred  region  of  higher  N  value  in  the  craniocaudal  (cc)  and 
oblique  (ob)  views  of  the  left  (1)  breast.  No  evident  structures  are  visible  in  the  N-images  for  the 
right  breast.  Figure  2  reports  the  N”  image.  The  cancer  appears  now  with  a  much  better  defined 
shape  in  the  left  breast,  while  new  structures  become  evident  in  the  right  and  left  breast 
mammograms.  In  particular,  we  assign  most  of  the  structures  that  appear  in  the  right  breast  N” 
images  to  blood  vessels. 

Task  l.d. 

During  year  01,  we  have  developed  a  perturbation  approach,  based  on  diffusion  theory,  to  assign  an 
absorption  perturbation  to  each  optical  inhomogeneity  detected  in  the  breast.  This  perturbation 
approach  is  based  on  the  fact  that  the  absorption  perturbation  (to  within  an  unknown 
wavelength  independent  factor  C)  can  be  written  as: 
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where  pj0'  is  the  average  reduced  scattering  coefficient  of  healthy  breast  tissue,  J0  is  the  unperturbed 
optical  intensity,  and  A I  is  the  measured  intensity  perturbation.  Finally,  we  have  fitted  the  spectrum 
of  the  absorption  perturbation  with  a  linear  combination  of  the  extinction  spectra  of  oxy-  and 
deoxy-hemoglobin.  By  considering  the  oxy-  and  deoxy-hemoglobin  concentrations  ([Hb02]  and 
[Hb],  respectively)  as  fitting  parameters,  we  can  translate  the  spectral  information  into  the  tumor 
oxygenation  (given  by  [Hb02]/([Hb]+[Hb02])). 

During  year  02,  we  have  continued  the  spectral  study  by  expanding  the  analysis  performed 
during  year  01  in  two  different  ways: 

1)  We  have  observed  that  in  the  perturbation  limit,  the  peak  value  of  the  second  derivative  of  N 
(N”max)  is  proportional  to  the  absorption  perturbation  introduced  by  the  tumor.  Therefore, 
we  have  performed  the  fit  of  the  spectrum  of  N”(A,)  (as  described  above)  to  find  an 
oxygenation  index  (01).  This  oxygenation  index,  which  in  a  real  perturbation  case  would 
coincide  with  the  hemoglobin  saturation  at  the  tumor,  is  only  qualitatively  related  to  the 
oxygenation  because  real  tumors  are  typically  beyond  the  validity  of  perturbation  theory. 
Of  course,  we  only  calculate  the  oxygenation  index  (01)  at  the  pixels  where  the  second 
derivative  is  negative  (gray  areas  in  Fig.  2).  If  we  now  identify  with  false  colors  the  regions 
with  lowest  values  of  oxygenation  index,  we  find  that  we  may  be  able  to  discriminate  blood 
vessels  from  cancer.  This  is  the  case  for  our  representative  patient,  for  whom  we  report  in 
Fig.  3  the  false  color  OI  images.  In  fact,  the  lowest  value  (blue  in  this  color  representation) 
only  occurs  at  the  cancer  location.  We  are  now  in  the  process  of  testing  this  new  idea  on  a 
subset  of  optical  mammogram,  to  extend  it  subsequently  to  the  whole  dataset. 

2)  We  have  found  that  the  right-hand  side  of  Eq.  (1)  is  not  simply  proportional  to  the 
absorption  perturbation  (ApJ  once  we  are  not  within  the  limits  of  perturbation  theory. 
However,  we  have  found  that  (A/)max//0  is  indeed  only  a  function  of  the  product  ps0'Apa  even 
in  cases  beyond  the  realm  of  perturbation  theory  (under  the  assumption  that  Ap/  =  0,  as 
approximately  true  for  breast  tumors).  We  have  found  this  result  by  calculating  -(A/)max//0 
for  the  case  of  a  large  sphere  (not  a  small  perturbation)  and  various  values  of  ps0'  and  Apa. 
The  fact  that  -(A7)max//0  only  depends  on  the  product  ps0'Apa  is  illustrated  in  Fig.  4(a),  while 
Fig.  4(b)  shows  the  independence  (to  a  good  approximation)  of  -(A/)max//0  on  the  background 
absorption  pa0.  This  has  prompted  the  development  of  a  new  idea,  which  consists  of 
identifying  two  wavelengths  Xj  and  X2  such  that  (A 7)max//0  is  the  same  at  the  two 
wavelengths.  It  turns  out  that  because  of  the  shape  of  the  hemoglobin  spectrum  and  the 
spectrum  of  typical  breast  tissue,  this  condition  can  be  fulfilled  to  a  good  approximation 
using  two  wavelengths  that  are  at  least  50  nm  apart  (otherwise  the  approach  has  a  large 
error).  The  basic  idea  is  that  at  these  two  wavelengths,  since  (A7)max//0  is  the  same,  then  also 
the  product  ps0'Apa  is  the  same,  and  therefore  we  can  write: 

AM*n)  =  ^0(^2)  (2) 

Apa(^2)  p^(Xj)  ' 

This  translates  a  measurement  of  the  spectral  ratio  of  the  background  reduced  scattering 
(which  can  be  performed  with  frequency-domain  methods)  into  a  measurement  of  the 
spectral  ratio  of  the  absorption  perturbation,  which  is  all  that  is  needed  for  a  quantitative 
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measurement  of  the  tumor  oxygenation.  Even  though  the  current  dataset  does  not  allow  us 
to  identify  the  two  wavelengths  X,  and  X2  from  a  highly  resolved  spectrum  (only  four 
wavelengths  are  available),  we  can  still  try  to  extrapolate  the  data  at  the  four  wavelengths  to 
a  continuous  spectrum  to  obtain  indications  on  the  applicability  of  our  idea  to  clinical  data. 
At  any  rate,  this  new  approach  is  in  line  with  the  broad  objective  of  this  project  of  guiding 
the  design  of  new  optical  instrumentation  for  breast  cancer  detection. 


Eao  (cm1) 


Fig.  4.  Theoretical  calculations  of  the  dependence  of  -(A/)max//0  for  a  spherical  inclusion  on  the  background  optical 
properties  and  on  the  absorption  perturbation.  The  scattering  perturbation  is  assumed  to  be  zero,  because  the  scattering 
contrast  of  tumors  is  typically  much  smaller  than  their  absorption  contrast. 
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Task  2 

We  have  performed  a  laboratory  study  that  was  divided  into  three  segments:  (1)  preparation  and 
optical  characterization  of  inclusions  to  simulate  tumors,  (2)  preparation  and  optical 
characterization  of  the  background  material  (a  highly-scattering,  liquid,  infinite  medium),  and  (3) 
experimental  procedure  to  investigate  that  A/7/0|  depends  solely  on  the  product  p's0Apa  (as 

theoretically  predicted  in  the  case  of  a  spherical  inclusion)  for  the  case  of  cylindrical  and 
irregularly-shaped  objects  embedded  in  a  uniform  turbid  medium.  Segments  (1)  and  (2)  pertain  to 
Task  2a,  while  segment  (3)  pertains  to  tasks  2b  and  2c. 

Task  2a 

(1)  Inclusions  were  comprised  of  a  mixture  of  General  Electric  silicones,  model  No. 
RTV615  (clear)  and  model  No.  RTV11  (white).  The  white  silicone  was  used  as  the  scattering 
material  and  black  India  ink  was  used  as  the  absorber.  Two  different  mixtures  were  prepared  with 
different  optical  properties,  (the  first  with  pa  ~  0.05-0.08  cm'1  and  p/  ~  9-12  cm'1  and  the  second 
with  pfl  ~  0.1 1-0.19  cm'1  and  p/  ~  9-1 1  cm'1  over  the  wavelength  range  of  752-840  nm  considered 
by  us)  to  cover  the  range  of  Apa  from  0.01  cm'1  to  0.2  cm'1  that  is  representative  of  breast  lesions  in 
vivo.  The  optical  properties  were  measured  using  frequency-domain  data.  The  inclusion  material 
was  formed  into  irregular  shapes  (by  cutting  cylindrical  shapes  with  razor  blades)  and  into  a  10  cm- 
long  cylinder  using  the  two  different  mixtures,  thus  creating  two  different  sets  of  the  same  sizes  and 
shapes.  The  cylinder  had  a  diameter  of  1.0  cm,  while  the  irregular  shapes  were  created  to  have  the 
same  volumes  as  spheres  with  1.0  cm  diameter  (V=  1.4  cm3)  and  1.4  cm  diameter  ( V  =  0.5  cm3), 
respectively. 

(2)  The  background  medium  was  comprised  of  1  liter  of  Liposyn  (10%  lipids  content)  to  8 
Liters  of  deionized  water  and  was  contained  inside  a  rectangular  vessel  with  dimensions  of  27  cm 
width,  32  cm  length,  and  13  cm  height.  Frequency-domain  measurements  of  pa0  over  the 
wavelengths  used  by  us  ranged  from  0.006  cm'1  to  0.04  cm'1  and  those  of  pi0'  ranged  from  9  cm'1  to 
15  cm1.  These  value  of  the  optical  properties  match  typical  near-infrared  absorption  and  scattering 
coefficients  in  breast  tissue  (pa(breast)  ~  0.03  cm'1  and  p/(breast)  ~  12  cm"1)  [1]. 

Tasks  2b  and  2c 

(3)  The  experimental  procedure  to  verify  that  A I  /  70|max  depends  on  the  product  p's0Apfl  was  based 

on  the  arrangement  shown  in  Fig.  5.  A  frequency-domain,  near-infrared  spectrometer  (ISS,  Inc., 
Champaign,  IL,  Model  No.  96208)  housed  and  controlled  the  laser  sources  and  optical  detectors 
used  to  obtain  the  measurements.  The  sources  were  laser  diodes  at  six  discrete  wavelengths  (752, 
778,  786,  813,  830  and  840  nm)  and  the  optical  detector  was  a  photomultiplier  tube  (Hamamatsu 
Photonics  R928).  The  sources  were  intensity-modulated  at  a  frequency  of  1 10  MHz,  multiplexed  at 
a  rate  of  27  Hz  for  discrimination  at  the  detector,  and  were  coupled  to  400  pm  core-diameter  optical 
fibers  that  were  collected  into  a  fiber  bundle  with  rectangular  cross-section  1 .2  mm  in  width  and 
3  mm  in  length.  The  detector  was  coupled  to  another  fiber  bundle  of  circular  cross-section,  3  mm 
internal  diameter.  The  ends  of  the  illumination  and  collection  fiber  bundles  were  fully  immersed  in 
the  Liposyn  mixture  to  simulate  an  infinite  medium.  The  source  fiber  and  the  collection  fiber  were 
then  moved  in  tandem,  remaining  collinear  all  the  time,  at  a  rate  of  0.65  mm/second.  Data  was 
acquired  at  every  0.5  seconds  providing  a  measurement  every  325  pm  over  a  scanning  range  of 


11 


14  cm.  The  source-detector  separation  for  these  experiments  remained  a  constant  6  cm,  which  is 
representative  of  the  thickness  of  a  slightly  compressed  breast.  Inclusions  of  different  sizes  and 
shapes  were  then  suspended  in  the  medium  equidistant  from  the  source  and  detector  using  Pasteur 
pipettes  to  hold  them  in  place.  The  pipettes  were  filled  with  background  medium  to  reduce  the 
optical  perturbation  they  might  cause. 


Fig.  5.  Experimental  setup  for  the  optical  measurements  in  a  tissue-like  strongly  scattering  medium  containing  tumor- 
simulating  inclusions.  MUX:  multiplexer,  PMT:  photo-multiplier  tube,  RF:  radio  frequency. 


Results 

The  experimental  results  for  A/  //0|max  as  a  function  of  the  product  pj0'Apa  are  reported  in  Fig.  6(a) 

for  the  irregularly-shaped  inclusions,  and  in  Fig.  6(b)  for  the  cylindrical  inclusion.  In  both  cases, 
the  fact  that  the  experimental  data  of  -A///0|mgx  for  a  range  of  values  of  pj0'  and  Apa  are 

distributed  along  a  single  curve  as  a  function  of  the  product  |aj0'Ap.fl  is  indicative  of  the  dependence 
of  A II  /0|max  on  the  product  pl0'Apa.  In  Figs.  6(a)  and  6(b),  the  functions  of  |ai0'A|aa  are  indicated 

by  arbitrary  continuous  lines.  The  small  deviations  of  the  experimental  points  from  the  arbitrary 
smooth  lines  are  assigned  to  experimental  errors  and  to  the  approximate  fulfillment  of  the  condition 
Ap/  =  0.  These  experiments  confirm  our  hypothesis  that  A/  / /0  |max  is  only  a  function  of  the 

product  mV  Aik  even  for  objects  with  a  relatively  large  size,  cylindrical  or  irregular  shape,  and 
relatively  high  absorption  contrast. 
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Fig.  6.  Experimental  results  of  the  dependence  of  -(A I)mJI0  on  the  product  of  the  background  scattering  coefficient  and 
the  absorption  perturbation,  (a)  Irregularly  shaped  inclusions,  (b)  cylindrical  inclusion. 


The  relevance  of  these  findings  is  illustrated  in  Fig.  7,  which  shows  that  we  can  accurately  measure 
the  oxygen  saturation  of  tumors  (modeled  here  by  spherical  inclusions)  even  if  their  size  and 
absorption  contrast  is  well  beyond  the  limits  of  perturbation  theory.  We  have  computed  the  values 
of  M /7olmax  at  nine  wavelengths  (680,  700,  720,  730,  758,  780,  800,  840,  and  880  nm)  for  a 

spherical  inclusion  in  a  case  that  mimics  a  realistic  condition  in  optical  mammography.  The 
background  scattering  and  absorption  spectra  were  set  equal  to  reported  spectra  for  healthy  breast 
tissue  [2],  we  set  the  hemoglobin  concentration  of  the  embedded  lesion  to  a  value  60  pM  higher 
than  that  in  the  background,  and  we  varied  the  hemoglobin  saturation  of  the  sphere  over  the  full 
range  0-100%.  We  have  compared  the  measurements  of  hemoglobin  saturation  obtained  (l)by 
using  perturbation  theory  [i.e.  Eq.  (1)]  at  two  fixed  wavelengths  and  (2)  by  using  our  method  [i.e. 
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Eq.  (2)  with  a  linear  correction  that  takes  into  account  deviations  from  the  condition 
(A 2)max//0(A-i)  =  at  two  wavelengths  that  depend  on  the  oxygenation  of  the  inclusion. 

Both  approaches  require  knowledge  of  the  background  reduced  scattering  coefficient,  which  can  be 
measured  with  frequency-domain  methods.  The  results  of  the  perturbation  analysis  for  a  number  of 
wavelength  pairs  are  shown  in  Fig.  7(a)  for  a  1 .4  cm-diameter  sphere  (on  the  mid-line  between 
source  and  detector),  in  Fig.  7(c)  for  a  3.0  cm-diameter  sphere  (on  the  mid-line  between  source  and 
detector),  and  in  Fig.  7(e)  for  a  1.4  cm  diameter  sphere  off  the  mid-line  by  1.5  cm.  The  results  of 
our  new  method  are  shown  in  Fig.  7(b)  (sphere  diameter:  1.4  cm;  on  the  mid-line),  Fig.  7(d)  (sphere 
diameter:  3.0  cm;  on  the  mid-line),  and  Fig.  7(f)  (sphere  diameter:  1.4  cm;  1.5  cm  off  the  mid-line). 
These  results  show  that  our  new  method  can  achieve  accurate  measurements  of  the  oxygenation  of 
spherical  regions  over  a  wide  range  of  oxygenation  values  (in  particular  from  40  to  100%)  and 
independent  of  the  size  and  location  of  the  sphere. 

Additional  research 

Some  of  our  results  on  the  measurements  of  oxygenation  levels  from  spectral  measurements  have 
been  developed  in  conjunction  with  an  effort  to  measure  the  oxygenation  of  the  venous 
compartment  in  the  skeletal  muscle  of  animal  models  (piglets)  and  human  subjects.  Our  approach 
and  findings  in  this  area  are  reported  in  the  manuscripts  enclosed  with  this  report. 


KEY  RESEARCH  ACCOMPLISHMENTS 

•  Development  of  a  new  scheme  of  image  processing  based  on  the  second  spatial  derivative  of 
the  N  images. 

•  Computation  of  a  tumor  oxygenation  index  (OI)  based  on  the  spectral  dependence  of  the 
second  derivative  of  N  at  the  tumor. 

•  Combination  of  the  above  two  points  into  a  composite  false-color  breast  image  containing 
structural  information  (from  the  second  derivative  of  N)  and  functional  information  (OI). 

•  Development  of  a  new  method  to  measure  the  tumor  oxygenation  from  an  appropriate 
choice  of  a  wavelength  pair  that  depends  on  the  tumor  oxygenation. 

•  Theoretical  tests  of  the  new  method  using  the  analytical  solution  for  a  spherical  inclusion 
(tumor)  embedded  in  a  turbid  medium  (breast  tissue). 

•  Experimental  tests  of  the  new  method  on  breast-simulating  turbid  media. 
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St02  (%)  -  fixed  X 


Calculations  using  a  fixed  wavelength  pair  over 
the  full  St02  range. 


Calculations  using  wavelength  pairs 
appropriately  chosen  for  each  St02  value. 


Fig.  7  (previous  page).  Analysis  of  the  effectiveness  of  our  new  method  for  measuring  the  tumor  oxygenation  (right 
panels)  with  respect  to  first  order  perturbation  theory  using  fixed  wavelength  pairs  (left  panels).  The  calculated 
oxygenation  values  (y  axes)  are  correlated  with  the  actual  oxygenation  values  ( x  axes)  used  in  the  theoretical  calculation 
of  the  data. 


REPORTABLE  OUTCOMES 

Manuscripts  (three  copies  of  these  manuscripts  are  enclosed  with  this  report) 

S.  Fantini  and  M.  A.  Franceschini,  “Frequency-Domain  Techniques  for  Tissue  Spectroscopy  and 
Imaging,”  in  Optical  Biomedical  Diagnostics,  V.  V.  Tuchin,  Ed.,  (SPIE  Press,  Bellingham,  WA, 
2001),  Chapter  7,  ( accepted  for  publication). 

M.  A.  Franceschini,  D.  A.  Boas,  A.  Zourabian,  S.  G.  Diamond,  S.  Nadgir,  D.  W.  Lin,  J.  B.  Moore, 
and  S.  Fantini,  “Near-Infrared  Spiroximetry:  Non-Invasive  Measurement  of  Venous  Saturation  in 
Piglets  and  Human  Subjects,”  J.  Appl.  Physiol.,  (accepted  for  publication). 

M.  A.  Franceschini,  A.  Zourabian,  J.  B.  Moore,  A.  Arora,  S.  Fantini,  and  D.  A.  Boas,  “Local 
Measurement  of  Venous  Saturation  in  Tissue  with  Non-Invasive,  Near-Infrared  Respiratory- 
Oximetry,”  Proc.  SPIE  4250,  1 64- 1 67  (200 1 ). 

Two  manuscripts  are  currently  in  preparation.  They  will  report  our  results  in  the  areas  of: 

1)  Quantitative  oximetry  of  breast  tumors  with  near-infrared  spectral  imaging; 

2)  Second  derivative  imaging  combined  with  oxygenation  imaging. 


CONCLUSIONS 

We  have  developed  a  new  image  processing  method  (second  derivative  imaging)  to  better  display 
the  optical  mammograms.  This  new  method  results  in  further  improvement  of  contrast  and 
detection  of  breast  tumors  with  respect  to  the  previous  N-imaging  approach.  We  have 
complemented  this  improved  structural  imaging  approach  with  a  functional  imaging  approach 
based  on  an  oxygenation  index  computed  from  spectral  features  of  the  second  derivative  images. 
We  have  also  refined  the  perturbative  approach  for  spectral  analysis  developed  during  year  01. 
This  refinement  has  led  to  a  new  method  for  measuring  the  tumor  oxygenation  with  an  excellent 
accuracy  (from  theoretical  data)  over  the  range  of  oxygenation  40-100%. 

The  relevance  of  our  results  is  that  single  wavelength  images  that  provide  structural 
information  can  be  superimposed  with  false  color  oxygenation  images  that  provide  functional 
information.  Our  new  approach  has  the  potential  of  quantifying  the  tumor  oxygenation,  so  that 
optical  discrimination  of  benign  and  malignant  tumors  may  become  feasible. 
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ABSTRACT 

We  present  preliminary  results  of  non-invasive,  near-infrared  measurements  of  venous  saturation  (Sv02)  on  the  leg  muscle  of 
three  anesthetized  piglets.  We  have  quantified  the  local  Sv02  by  analyzing  the  optical  spectrum  of  the  amplitude  of  the 
absorption  oscillations  synchronous  with  breathing.  To  induce  a  variation  in  the  muscle  Sv02,  we  performed  measurements 
during  a  protocol  involving  a  cyclic  change  in  the  fraction  of  oxygen  inspired  by  the  piglet  over  the  range  10-100%  (by 
volume).  In  all  three  piglets,  we  have  found  a  good  agreement  between  the  Sv02  values  measured  non-invasively  with  near- 
infrared  spectroscopy  (NIRS)  and  those  measured  invasively  by  the  analysis  of  venous  blood  samples. 

Keywords:  Near-infrared,  tissue  spectroscopy,  frequency-domain,  pulse  oximetry,  venous  saturation. 


1.  INTRODUCTION 

The  hemoglobin  concentration  in  blood-perfused  tissues  shows  temporal  fluctuations  that  are  induced  by  physiological  and 
metabolic  activity.  For  example,  the  arterial  pulsation  associated  with  the  systolic-diastolic  blood  pressure  variation  induces 
hemoglobin  concentration  oscillations  at  the  heartbeat  frequency.  These  oscillations  are  exploited  by  pulse  oximeters  to 
measure  the  arterial  saturation.1-2  Other  fluctuations  on  a  longer  time  scale  of  seconds  or  tens  of  seconds  are  associated  with 
vasomotion  activity,  heart  rate  fluctuations,  and  breathing.3'6  Near-infrared  tissue  spectroscopy  (NIRS),  being  exquisitely 
sensitive  to  the  absorption  of  hemoglobin,  is  an  effective  tool  to  measure  the  hemoglobin  concentration  in  tissues  and  its 
temporal  fluctuations.  Furthermore,  one  can  optically  discriminate  the  oxy-  and  deoxy-hemoglobin  species  because  of  their 
different  near-infrared  absorption  spectra. 

In  this  article,  we  present  an  optical  study  of  the  hemoglobin  fluctuations  at  the  breathing  frequency  in  the  leg  muscle  of 
newborn  piglets.  Under  the  hypothesis  that  the  venous  compartment  gives  the  dominant  contribution  to  the  hemoglobin 
oscillations  at  the  breathing  frequency,  we  have  used  this  oscillatory  component  to  measure  the  venous  saturation.  This 
hypothesis  was  originally  reported  by  Wolf  et  al?  This  measurement  of  venous  saturation  is  local,  non-invasive,  can  be  done 
continuously,  and  does  not  require  any  perturbation  such  as  cuff  occlusion,  dye  injection,  or  gas  inhalation. 
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2.  METHODS 


2.1.  Tissue  spectrometer 

The  optical  measurements  were  performed  with  a  multi-channel  tissue  spectrometer  operating  in  the  frequency-domain 
(Model  96208,  ISS,  Inc.,  Champaign,  IL).8'9  For  the  measurements  reported  here,  we  employed  two  photomultiplier  tube 
(PMT)  detectors  and  eight  laser  diodes  emitting  at  633,  670,  691,  752,  778,  786,  830,  and  840  nm.  The  eight  light  sources, 
which  were  intensity  modulated  at  a  frequency  of  1 10  MHz,  were  electronically  multiplexed  at  a  rate  of  100  Hz  to  time-share 
the  two  parallel  PMT  detectors.  As  a  result,  a  complete  acquisition  cycle  over  the  eight  wavelengths  took  80  ms.  The  laser 
diodes  and  the  PMT's  were  all  coupled  to  fiber  optics.  The  eight  individual  source  fibers,  each  400  pm  in  internal  diameter, 
were  closely  arranged  into  a  fiber  bundle  having  a  rectangular  cross-section  of  3.5  x  2.0  mm2.  The  detector  fiber  bundles 
were  1 .5  mm  in  internal  diameter  (we  also  tested  3.0  mm  fiber  bundles  and  we  found  no  significant  difference  with  respect  to 
the  results  obtained  with  the  1.5  mm  bundles).  The  optical  fibers  were  placed  in  contact  with  the  skin  of  the  piglet  by  means 
of  a  flexible,  plastic  probe.  This  probe  arranged  the  tips  of  the  source  and  detector  fiber  bundles  along  a  line  with  the  two 
detector  bundles  at  distances  of  1 .0  and  2.0  cm,  respectively,  from  the  source  bundle.  The  data  processing  and  instrument 
control  were  performed  by  software  running  under  Windows  98.  This  software  allowed  for  the  real-time  display  of  the 
average  intensity,  amplitude,  and  phase  data. 

2.2.  Measurement  protocol  1 

We  measured  three  piglets.  Each  piglet  was  anesthetized  (3-4%  Isoflurane  inhalation)  and  intubated.  The  animal  was  not 
ventilated  and  breathed  freely  throughout  the  experiment.  Femoral  cut-downs  were  performed  into  the  left  inferior  femoral 
vein  and  artery  to  insert  catheters  for  periodic  arterial  and  femoral  venous  blood  gas  monitoring.  The  optical  probe  was 
always  located  on  the  right  (non-catheterized)  posterior  leg,  with  one  exception  (piglet  #3)  where  data  were  also  collected  on 
the  catheterized  leg  (left)  for  comparison.  On  the  first  piglet  (14  days  old,  3.75  kg  body  weight)  we  performed  two 
oxygenation  cycles.  Each  cycle  consisted  of  varying  the  fraction  of  inspired  oxygen  (Fi02)  approximately  every  four  minutes 
through  the  values  of  40%,  15%,  10%,  and  100%  (by  volume).  During  each  specific  value  of  Fi02  we  acquired  about  3,000 
optical  data  points  (4  min  x  [60  s/min]  /  [80  ms/data  point]).  On  the  second  piglet  (16  days  old,  5.5  kg  body  weight)  we 
repeated  the  Fi02  cycle  four  times  (after  the  first  two  cycles  we  moved  the  optical  probe  to  a  different  location  on  the  same 
leg).  The  only  difference  with  respect  to  the  Fi02  cycle  performed  on  piglet  #1  was  that  on  piglet  #2  we  maintained  each 
particular  Fi02  level  for  about  six  minutes  (i.e.  about  4,500  data  points).  On  the  third  piglet  (15  days  old,  5.39  kg  body 
weight)  we  performed  three  Fi02  cycles  where  the  02  volume  fraction  assumed  the  sequential  values  of  40%,  20%,  17.5%, 
15%,  12.5%  10%,  and  100%.  Each  Fi02  value  was  maintained  for  about  six  minutes.  In  all  three  piglets,  an  invasive 
measurement  of  Sv02  from  a  femoral  vein  blood  sample  was  performed  at  the  end  of  each  Fi02  interval  using  a  conventional 
blood  gas  analyzer  (Instrumentation  Laboratory,  model  13040-11  pH/Blood  Gas  Analyzer).  Throughout  the  experiment,  we 
monitored  the  arterial  pulsation,  heart  rate  and  arterial  saturation  (Sa02)  with  a  commercial  pulse  oximeter  (Nellcor,  N-200) 
attached  to  the  piglet's  right  foot,  and  the  respiratory  motion  with  a  strain  gauge  belt  (Sleepmate/Newlife  Technologies, 
Resp-EZ)  placed  around  the  piglet's  chest.  The  readings  of  the  pulse  oximeter  and  the  strain  gauge  were  recorded 
continuously  and  simultaneously  with  the  collection  of  near-infrared  data.  Motion  artifacts  in  the  optical  data  were 
minimized  by  securing  the  piglet's  leg  to  the  operating  table. 

2.3.  Data  analysis 

We  used  the  modified  Beer-Lambert  law  approach10  to  translate  the  temporal  intensity  changes  into  variations  in  the  tissue 
absorption,  under  the  assumption  that  the  tissue  scattering  is  slowly-varying.  We  continuously  measured  the  optical 
pathlength  using  our  frequency-domain  data  and  diffusion  theory.  The  mean  pathlength  was  measured  at  the  end  of  each 
Fi02  period  by  quantifying  the  average  (over  1,000  points,  or  80  s)  absorption  and  reduced  scattering  coefficients  of  the 
tissue  using  the  frequency-domain  multi-distance  method.11  In  this  case,  the  multi-distance  scheme  was  realized  by  the  two 
detector  fiber  bundles  at  two  different  distances  (1.0  and  2.0  cm)  from  the  single  source  fiber  bundle.  The  different  fiber 
throughput  and  detector  sensitivity  for  the  two  detector  channels  was  accounted  for  by  a  preliminary  calibration  measurement 
on  an  optically  characterized  synthetic  tissue-like  sample.  The  validity  of  the  calibration  procedure  was  verified  at  the  end  of 
each  measurement  session  by  re-positioning  the  optical  probe  on  the  calibration  sample.  The  measured  values  of  the  optical 
pathlengths  in  the  piglet's  leg  at  the  eight  wavelengths  were  updated  every  time  the  Fi02  was  changed.  This  long  integration 
time  for  the  pathlength  measurements  (80  s)  realized  a  low-pass  filter  that  minimized  the  time-varying  contributions  from  the 
hemoglobin  oscillations  resulting  from  the  arterial  pulsation  and  breathing.  Furthermore,  the  multi-distance  measurement 
scheme  for  the  pathlength  measurement  provided  some  level  of  spatial  averaging  as  well.  By  contrast,  the  optical  data  for  the 
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measurement  of  venous  saturation  were  acquired  with  an  80  ms  temporal  resolution  and  using  a  single  source-detector 
distance  ( 1  cm). 

To  compute  the  venous  saturation,  we  measured  the  amplitude  of  the  absorption  oscillations  at  the  breathing  frequency  at  the 
eight  wavelengths  considered.  We  quantified  the  amplitude  of  the  breathing-induced  absorption  oscillations  by  taking  the 
sum  of  the  amplitudes  of  the  fast  Fourier  transform  of  Aun  over  the  breathing  band  or  its  second  harmonic.  We  evaluated  the 
fast  Fourier  transform  of  A(ifl  over  256  data  points,  corresponding  to  a  time  trace  20.5  s  long,  to  achieve  reliable  and 
reproducible  spectra  from  a  number  of  breathing  periods.  Furthermore,  we  averaged  about  800  successive  Fourier  transforms 
(each  computed  from  a  data  set  shifted  by  one  data  point  with  respect  to  the  previous  one),  so  that  the  total  number  of  data 
points  resulting  in  a  single  venous  saturation  reading  was  on  the  order  of  1,000  (corresponding  to  a  train  of  data  80  s  long). 
This  80  s-long  data  set  was  chosen  to  be  at  the  end  of  each  Fi02  period.  The  spectrum  of  the  amplitude  of  the  absorption 
oscillations  synchronous  with  breathing  was  fitted  to  a  linear  combination  of  the  oxy-  and  deoxy-hemoglobin  extinction 
spectra.  The  fitted  parameters  were  the  oscillatory  concentrations  of  oxy-hemoglobin  ([Hb02]osc)  and  deoxy-hemoglobin 
([Hb]osc).  The  hemoglobin  saturation  of  the  hemoglobin  compartment  oscillating  synchronously  with  breathing  is  then  given 
by  [HbCbLsc  /([Hb02]osc  +[Hb]0SC).  We  observe  that  this  procedure  can  be  applied  to  an  oscillatory  hemoglobin  component  at 
any  frequency.  At  the  heartbeat  frequency,  this  procedure  yields  the  arterial  saturation.12  Here,  we  apply  it  to  the  breathing 
frequency  to  obtain  the  venous  saturation. 


3.  RESULTS 

Figure  1  shows  typical  temporal  traces  (panel  (a))  and  amplitude  spectra  (panel  (b))  of  the  relative  oxy-  and  deoxy- 
hemoglobin  concentrations  measured  on  the  piglet's  leg.  The  temporal  traces  of  [Hb]  and  [Hb02]  are  obtained  by  fitting  the 
measured  absorption  spectrum  at  each  time  with  a  linear  combination  of  the  oxy-  and  deoxy-hemoglobin  extinction  spectra. 
Figure  1(b)  shows  the  values  of  the  venous  saturation  measured  from  the  fundamental  (0.63  Hz;  Sv02  =  78%)  and  from  the 
second  harmonic  (1.26  Hz;  Sv02  =  75%)  of  the  breathing  period.  Figure  1(b)  also  shows  the  [Hb02]  and  [Hb]  peaks  at  the 
heartbeat  frequency  (2.54  Hz)  from  which  we  have  been  able  to  measure  the  arterial  saturation  (Sa02  =  99%).  By 
comparison,  the  invasive  measurement  of  venous  and  arterial  saturation  with  the  blood  gas  analyzer  gave  values  of  77.5% 
and  99%,  respectively.  The  oscillatory  components  of  the  [Hb]  and  [Hb02]  traces  at  the  breathing  and  heartbeat  frequencies 
are  readily  identifiable  in  Fig.  1(a),  which  also  shows  the  breathing  wave  recorded  by  the  strain  gauge  placed  around  the 
piglet's  chest.  In  a  number  of  cases,  we  found  more  reliable  results  for  Sv02  by  analyzing  the  absorption  peaks  at  the  second 
harmonic  rather  than  at  the  fundamental  of  the  breathing  rate.  In  Fig.  1  we  show  the  time  traces  and  the  amplitude  spectra  of 
the  oxy-  and  deoxy-hemoglobin  concentrations  because  they  are  easier  to  visualize  than  the  eight  absorption  traces  and 
spectra.  We  note  that  computing  the  venous  saturation  from  the  oxy-  and  deoxy-hemoglobin  traces  is  equivalent  to  using  the 
eight  absorption  traces,  provided  that  the  absorption  oscillations  at  the  eight  wavelengths  are  in  phase.  This  follows  from  the 
fact  that  both  the  Fourier  transform  and  the  fitting  procedure  to  translate  the  absorption  changes  into  hemoglobin 
concentration  changes  are  linear  operations. 

Figure  2  shows  the  comparison  between  the  venous  saturation  values  measured  non-invasively  with  near-infrared 
spectroscopy  (open  squares)  and  invasively  with  a  gas  analysis  of  a  venous  blood  sample  (filled  circles).  In  Fig.  2,  we  report 
the  results  for  the  first  two  Fi02  cycles  in  piglet  #1  (panel  (a))  and  piglet  #2  (panel  (b)),  and  for  the  last  Fi02  cycle  in  piglet 
#3.  Figure  2  also  shows  the  temporal  evolution  of  Fi02  during  the  measurement. 

Figure  3  summarizes  our  results  on  the  three  piglets  examined.  The  shape  of  the  symbols  indicate  the  piglet  (circles:  piglet 
#1;  squares:  piglet  #2;  triangles:  piglet  #3),  while  the  filling  indicates  the  side  of  the  near-infrared  measurement  (filled:  right 
leg;  open:  left  leg). 


166 


Proc.  SP1E  Vol.  4250 


breathing  wave 


V. 

s 

-0.1 

3 

-O 

r—i 

C-4 

O 

-0.3 

X 

X 

U_ J 

-0.4 

o 

> 

-0.5 

•*— * 

-0.6 

•— < 

-0.7 

-0.8 

/^vAS\r/iAf‘ 


210 


212 


214  216 

time  (sec) 


218 


220 


s 

=L 

W' 

X> 

x 

o> 

> 

2 


(a) 


breathing  rate  arterial  pulse 


frequency  (Hz) 


Fig.  1.  (a)  Temporal  traces  and  (b)  amplitude  spectra  of  the  oxy-  and  deoxy-hemoglobin  concentrations  ([HbCT]  and 
[Hb],  respectively)  measured  in  piglet  #1  while  Fi02  was  set  to  40%.  Panel  (a)  also  shows  the  signal  recorded  by  the 
strain  gauge  belt  to  monitor  the  breathing  wave.  The  venous  saturation  (SvCF)  measured  from  the  amplitudes  of  the 
Fourier  transforms  at  the  fundamental  and  second  harmonic  of  the  breathing  frequency  are  shown  in  panel  (b),  which 
also  shows  the  arterial  saturation  (SaCK)  measured  from  the  peak  at  the  arterial  pulsation  frequency.  For  comparison, 
the  values  of  Sv02  and  Sa02  measured  invasively  by  the  blood  gas  analyzer  were  77.5%  and  99%,.  respectively. 
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piglet#! 


piglet  #2 


piglet  #3 


□  Sv02  from  NIRS  (non-invasive)  •  Sv02  from  blood  gas  analyzer  (invasive) 


Fig.  2.  Comparison  between  the  venous  saturation  values  measured  non-invasively  with  NIRS  (squares)  and  invasively 
from  blood  samples  (circles)  on  (a)  piglet  #1,  (b)  piglet  #2,  and  (c)  piglet  #3.  The  temporal  evolution  of  the  FiCF  level 
during  the  measurement  is  also  shown.  The  lines  that  connect  the  points  are  added  to  guide  the  eye.  The  error  bars 
indicate  standard  deviation  errors  from  the  approximately  800  Fourier  transforms  averaged  for  each  NIRS  reading  of 
SvOj. 
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i  Sv02  from  blood  gas  analyzer  (invasive)  (%) 

Fig.  3.  Overall  correlation  between  the  venous  saturation  values  measured  non-invasively  with  NIRS  (y-axis)  and 
invasively  by  a  blood  gas  analyzer  (x-axis).  The  circles  refer  to  piglet  #1,  the  squares  to  piglet  #2,  and  the  triangles  to 
piglet  #3.  The  open  triangles  indicate  two  FiCF  cycles  where  the  non-invasive  optical  data  were  acquired  on  the  same 
leg  (left)  that  was  catheterized  for  the  collection  .of  blood  samples  for  the  invasive  measurement  of  SvO> 


4.  DISCUSSION 

We  found  a  good  agreement  between  the  venous  saturation  values  measured  non-invasively  with  near-infrared  Fourier- 
oximetry  and  invasively  by  gas  analysis  of  venous  blood  samples.  This  can  be  seen  in  the  representative  traces  shown  in 
Fig.  2  as  well  as  in  the  overall  summary  shown  in  Fig.  3.  All  the  results  reported  here  are  based  on  the  analysis  of  the  peaks 
at  the  second  harmonic  of  the  breathing  frequency.  While  during  normoxia  and  hyperoxia  the  results  obtained  by  using  the 
fundamental  frequency  are  comparable  to  the  results  based  on  the  second  harmonic,  in  some  hypoxic  cases  the  calculations  at 
the  fundamental  frequency  yielded  negative  values  of  Sv02.  The  invasive  measurements  were  always  performed  on  the  left 
leg,  while  the  non-invasive  measurements  were  usually  performed  on  the  right  leg.  This  raises  a  question  as  to  whether  it 
makes  sense  to  compare  measurements  on  different  legs.  We  have  found  indications  that  under  our  experimental  conditions 
this  comparison  is  meaningful.  In  fact,  several  near-infrared  measurements  performed  on  the  same  leg  used  for  the  collection 
of  blood  samples  (left)  have  reproduced  the  measurements  on  the  right  leg  (see,  for  instance  the  open  symbols  in  Fig.  3). 
Furthermore,  we  did  not  find  significant  differences  by  moving  the  probe  to  different  locations  of  the  same  leg. 


5.  CONCLUSION 

We  have  reported  our  initial  results  of  venous  saturation  measurements  based  on  the  analysis  of  the  near-infrared  oscillations 
with  breathing.  The  non-invasive  assessment  of  venous  saturation  by  near-infrared  tissue  spectroscopy  complements 
previously  demonstrated  optical  measurements  of  tissue13'16  and  arterial12,17  saturation.  Therefore,  these  preliminary  results 
may  lead  to  the  design  of  a  single  instrument  capable  of  providing  simultaneous  non-invasive  readings  of  local  arterial, 
tissue,  and  venous  saturation. 
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[  David  W.  Lin,  John  B.  Moore,  and  Sergio  Fantini.  Near- 
infrared  spiroximetry:  noninvasive  measurements  of  venous 
;  saturation  in  piglets  and  human  subjects.  J  Appl  Physiol  91: 

;  000-000,  2001. — We  present  a  noninvasive  method  to  mea- 
'  ’  sure  the  venous .  oxygen  saturation  (|>v^  in  tissues  using 

near-infrared  spectroscopy  (NIRS).  This  method  is  based  on 
the  respiration-induced  oscillations  of  the  near-infrared  ab- 
AQ:2  sorption  in  tissues,  and  we  call  it  spiroximetry  (the  prefix 
yaspiro  means  respiration).  We  have  tested  this  method  in 
three  piglets  (hind  leg)  and  in  eight  human  subjects  (vastus 
.  medialis  and  vastus  lateralis  muscles).  In  the  piglet  study, 
}  we  compared  our  NIRS  measurements  of  the  Svg,  (SvA- 
NIRSreap)  with  the  Svg  of  blood  samples.  Sv^-NlRSreap  and 
Sv<|  of  blood  samples'" agreed  well  over  the  whole  range  of 
Svq2  considered  (20-95%).  The  two  measurements  showed 
an  average  difference  of  1.0%  and  a  standard  deviation  of  the 
f  difference  of  5.8%.  In  the  human  study,  we  found  a  good 
agreement  between  Svoj-NIRS^p  and  the  Svo2  values  mea¬ 
sured  with  the  NIRS  venous  occlusion  method.  Finally,  in  a 
'  "  preliminary  test  involving  muscle  exercise,  Svo2-NlRSresp 
showed  an  expected  postexercise  decrease  from  the  initial 
baseline  value  and  a  subsequent  recovery  to  baseline. 

tissue  spectroscopy;  frequency-domain;  pulse  oximetry;  he¬ 
ll  moglobin  saturation 
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the  possibility  of  using  light  to  measure  the  oxygen 
saturation  of  hemoglobin  in  vivo  has  been  explored 
since  the  1940s  (37).  The  feasibility  of  optical  blood 
oximetry  stems  from  the  oxygenation  dependence  of 
the  optical  spectrum  of  hemoglobin.  This  is  illustrated 
in  Fig.  1,  which  shows  the  absorption  spectra  of  100 
p-M  hemoglobin  for  oxygen  saturation  values  of  0,  20, 
40,  60,  80,  and  100%.  The  spectra  of  Fig.  1  were 


calculated  from  published  values  of  the  molar  extinc¬ 
tion  coefficients  of  oxyhemoglobin  (HbC>2)  and  deoxyhe¬ 
moglobin  (Hb)  (43,  53). 

Oxygen  saturation  of  the  pulmonary  capillary  blood 
in  rabbits  has  been  measured  by  using  dynamic  inva¬ 
sive  techniques  (48).  Near-infrared  light  in  the  wave¬ 
length  range  from  700  to  900  nm  results  in  a  sufficient 
penetration  depth  for  the  noninvasive  optical  monitor¬ 
ing  of  skeletal  muscle,  cerebral  gray  matter,  and  breast 
tissue.  As  a  result,  near-infrared  techniques  allow  a 
noninvasive  assessment  of  hemoglobin  saturation  for  a 
wide  range  of  applications,  such  as  the  study  of  muscle 
metabolism  (7,  9,  12,  29,  45),  the  diagnosis  of  vascular 
disorders  (2,  20,  32,  33,  44,  49),  functional  brain  imag¬ 
ing  (3,  10,  24,  30,  35,  50),  and  breast  cancer  detection 
(23,  28,  40,  42,  46). 

If  near-infrared  light  is  highly  sensitive  to  the  oxy¬ 
gen  saturation  of  hemoglobin,  then  its  large  penetra¬ 
tion  depth  inside  tissues  implies  that  the  arterial, 
venous,  and- capillary  compartments  all  contribute  to 
the  optical  signal.  The  average  hemoglobin  oxygen¬ 
ation  measured  with  near-infrared  spectroscopy 
(NIRS)  (19,  34,  41)  is  usually  referred  to  as  tissue 
oxygen  saturation  (StQ2).  St$2  values  are  assumed  to  be 
in  between  arterial  and  local  venous  saturation  values 
(Sa®2  and  Sv$,  respectively).  A  number  of  research 
studies  have  investigated  the  relationship  between  the 
near-infrared  (noninvasive)  measurement  of  Sto2  and 
the  values  of  Sao2  and  local  Svo2  measured  invasively 
from  drawn  blood  samples  (31,  51).  The  contribution  of 
the  arterial  compartment  to  the  noninvasive  optical 
signal  can  be  isolated  because  of  its  unique  temporal 
dynamics  associated  with  the  systolic-diastolic  blood 
pressure  variation  at  the  heartbeat  frequency  (1).  The 
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wavelength  (nm) 

Fig.  1.  Near-infrared  absorption  spectra  of  100  (iM  hemoglobin  con¬ 
centration  ([Hb-T|,  where-T,  stands  for  total)  for  different  values  of 
the  oxygen  saturation in  the  range  of  0-100%.  The  curve  for 
S02  =  0%  corresponds''!*)  the  deoxyhemoglobin  (Hb)  absorption 
spectrum,  whereas  the  curve  for  S02  =  100%  corresponds  to  the 
oxyhemoglobin  (Hb02)  absorption  spectrum.  These  spectra  have 
been  computed  from  published  spectra  of  the  molar  extinction  coef¬ 
ficients  of  Hb02  and  Hb  (43,  53). 

pulsatile  component  of  the  optical  signals  at  two  or 
more  wavelengths  at  the  heartbeat  frequency  is  used 
by  conventional  (I,  36)  or  self-calibrated  (21)  pulse 
oximeters  to  measure  the  Sa©2.  Sa8>2  is  a  parameter 
that  provides  information  about  the  ventilation  and 
the  oxygen  exchange  in  the  lungs.  In  contrast,  Svo2  is  a 
parameter  that  reflects  the  local  balance  between  blood 
flow  and  oxygen  consumption.  The  noninvasive  optical 
measurement  of  Svo2  is  complicated  by  the  fact  that 
the  isolation  of  the  contribution  of  the  venous  compart¬ 
ment  to  the  noninvasive  optical  signal  is  not  straight¬ 
forward.  There  are  no  clinical  devices  presently  capa¬ 
ble  of  monitoring  Svo2  noninvasively. 

A  number  of  experimental  approaches  have  been 
proposed  to  measure  Sv02  from  induced  local  changes 
in  the  venous  blood  volume.  For  instance,  proposed 
approaches  involve  a  venous  occlusion  in  a  limb  (13,  39, 
55,  56),  tilting  the  patient’s  head  down  by  15  degrees 
(47),  a  partial  jugular  vein  occlusion  (15,  54),  or  me¬ 
chanical  ventilation  (52).  In  all  these  approaches,  Svo2 
is  optically  measured  as  the  ratio  between  the  in¬ 
creases  in  the  Hb02  concentration  ([Hb02])  and  the 
total  hemoglobin  concentration  (equal  to  [Hb02]  + 
[Hb],  where  [Hb]  is  deoxyhemoglobin  concentration) 
induced  by  the  local  increase  in  venous  blood  volume. 
To  overcome  the  limitations  of  these  methods,  which 
can  either  be  applied  only  to  the  limbs  (venous  occlu¬ 
sion  method)  or  require  an  external  perturbation  (par¬ 
tial  jugular  vein  occlusion,  mechanical  ventilation,  and 
tilting  methods),  we  propose  an  alternative  approach 
that  is  an  extension  of  the  method  of  Wolf  et  al.  (52). 
This  approach  involves  no  external  perturbations  and 
is  applicable  to  subjects  who  are  breathing  either  spon¬ 
taneously  or  synchronously  with  a  metronome  set  at 
their  average  respiratory  frequency.  Furthermore,  this 
method  can  provide  continuous  and  real-time  monitor¬ 
ing  of  Svo2.  The  basic  idea  is  to  measure  Svo2  from  the 


amplitude  of  the  optically  measured  [HbCy  and  [Hb] 
oscillations  at  the  respiratory  frequency.  The  basic 
hypothesis,  originally  formulated  in  this  context  and 
tested  on  the  brain  of  mechanically  ventilated  infants 
by  Wolf  et  al.,  is  that  the  oscillatory  components  of 
[Hb02]  and  [Hb]  at  the  breathing  rate  are  mostly 
representative  of  the  venous  compartment.  Because 
the  venous  compliance  is  ~20  times  as  large  as  the 
arterial  compliance  (4),  a  given  change  in  the  blood 
pressure  in  the  veins  causes  a  venous  volume  change 
~20  times  as  large  as  the  arterial  volume  change 
corresponding  to  the  same  pressure  change  in  the  ar¬ 
teries.  During  normal  breathing,  the  inspiration  phase 
involves  a  decrease  in  the  intrathoracic  pressure  and 
an  increased  pressure  gradient  between  the  peripheral 
venous  system  and  the  intrathoracic  veins.  This  causes 
blood  to  be  drawn  from  the  extrathoracic  veins  into  the 
intrathoracic  vessels  and  heart  (26).  Because  of  the 
vein  valves,  venous  return  is  increased  more  by  inspi¬ 
ration  than  it  is  decreased  by  expiration  (38).  The  net 
effect  is  the  so-called  respiratory  pump  that  facilitates 
the  venous  return  from  the  periphery  by  the  respira¬ 
tion-induced  periodic  fluctuations  in  the  central  venous 
pressure  (38).  As  a  result  of  the  respiratory  pump,  the 
peripheral  venous  blood  volume  oscillates  at  the  respi¬ 
ratory  frequency,  decreasing  during  inspiration  and 
increasing  during  expiration. 

It  is  on  this  oscillatory  component  at  the  respiratory 
frequency  that  we  base  our  near-infrared  measure¬ 
ment  of  the  Svo2.  We  coin  the  term  spiroximeter  to 
indicate  an  instrument  for  measuring  the  Svo2  from 
respiration-induced  oscillations  in  the  venous  blood 
pressure  and  in  the  venous  volume  fraction  in  tissues. 
It  must  be  observed  that  respiration  may  also  induce 
perturbations  to  the  heart  rate  (respiratory  sinus  ar¬ 
rhythmia)  and  consequently  to  the  cardiac  output  and 
arterial  blood  pressure.  As  a  result,  the  arterial  com¬ 
partment  volume  may,  in  general,  also  oscillate  at  the 
respiratory  frequency;  thus  near-infrared  spiroximetry 
data  must  be  carefully  examined  to  guarantee  a  reli¬ 
able  reading  of  Svo2. 

We  report  a  validation  study  conducted  on  the  hind 
leg  of  three  piglets,  in  which  we  compared  the  near- 
infrared  measurements  of  Svo2  (Svo2-NIRS)  with  the 
Svo2  values  obtained  by  the  gas  analysis  of  venous 
blood  samples  (Svo2-blood).  To  show  the  applicability  of 
spiroximetry  to  human  subjects,  we  also  conducted  a 
preliminary  test  on  the  vastus  medialis  and  vastus 
lateralis  muscles  of  healthy  volunteers  at  rest  and 
postexercise. 

MATERIALS  AND  METHODS 

Tissue  spectrometer.  The  near-infrared  measurements 
were  performed  with  a  frequency-domain  tissue  spectrome¬ 
ter  (model  96208,  ISS,  Champaign,  IL)  (18,  25).  This  instru¬ 
ment  uses  two  parallel  photomultiplier  tube  detectors  that 
are  time  shared  by  eight  multiplexed  laser  diodes  emitting  at 
636,  675,  691,  752,  780,  788,  830,  and  840  nm,  respectively. 
The  frequency  of  intensity  modulation  is  110  MHz,  and 
heterodyne  detection  is  performed  with  a  cross-correlation 
frequency  of  5  kHz.  The  multiplexing  rate,  i.e.,  the  frequency 
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of  sequential  laser  switching,  is  100  Hz.  As  a  result,  50 
cross-correlation  periods  are  acquired  during  the  on  time  of 
each  laser  diode,  and  a  complete  acquisition  cycle  over  the 
eight  wavelengths  is  completed  every  80  ms.  The  laser  diodes 
and  the  photomultiplier  tubes  are  all  coupled  to  fiber  optics. 
The  eight  individual  illumination  fibers,  each  400  pm  in 
internal  diameter,  are  arranged  into  a  fiber  bundle  having  a 
rectangular  cross-section  of  3.5  X  2.0  mm2.  The  collecting 
circular  fiber  bundles  are  3.0  mm  in  internal  diameter.  The 
optical  fibers  are  placed  in  contact  with  the  skin  by  means  of 
a  flexible  plastic  probe.  The  optical  probe  arranges  the  tips  of 
the  illuminating  and  collecting  fiber  bundles  along  a  line, 
with  the  two  collecting  fiber  bundles  at  distances  of  1.0  and 
2.0  cm  from  the  single  illuminating  bundle.  In  some  cases,  we 
have  used  a  second  tissue  spectrometer  to  perform  simulta¬ 
neous  measurements  on  both  legs  (piglets  2  and  3)  or  at 
different  locations  on  the  same  leg  (human  subjects).  In  the 
second  tissue  spectrometer  (which  used  the  optical  probes  PL 
and  HVL  defined  below),  the  840-nm  laser  diode  was  replaced 
by  a  laser  diode  emitting  at  814  nm. 

Measurements  on  piglets.  We  performed  measurements  on 
three  piglets  that  were  15  ±  1  days  old  and  weighed  5  ±  1  kg. 
The  experimental  arrangement  for  the  piglet  measurements 
is  schematically  illustrated  in  Fig.  2.  The  piglets  were  anes¬ 
thetized  by  inhalation  of  3—4%  isoflurane  administered  by 
means  of  a  breathing  mask  applied  to  the  piglet’s  snout.  The 
animals  were  not  mechanically  ventilated,  and  they  breathed 
freely  throughout  the  experiment.  A  strain-gauge  belt  (Sleep- 
mate/Newlife  Technologies,  Resp-EZ)  was  placed  around  the 
piglet’s  thorax  to  continuously  monitor  the  respiratory  excur¬ 
sion.  A  pulse  oximeter  (Nellcor,  N-200)  continuously  recorded 
the  heart  rate  at  the  foot  of  the  right  hind  leg.  The  analog 
outputs  from  the  strain  gauge  and  the  pulse  oximeter  were 
fed  to  the  auxiliary  input  ports  of  the  tissue  spectrometer  for 
continuous  coregistration  of  optical  and  physiological  data.  A 
femoral  cutdown  was  performed  into  the  left  inferior  femoral 
vein  to  insert  a  catheter  for  periodic  blood  sampling.  The 
femoral  venous  blood  samples  were  run  through  a  commer¬ 
cial  blood-gas  analyzer  (Instrumentation  Laboratory,  model 


strain  gauge  belt 
(respiratory  excursion) 


1304  pH/blood-gas  analyzer)  to  obtain  invasive  readings  of 
Svo2-blood.  One  optical  probe  (identified  as  probe  PR)  was 
always  located  on  the  right  (noncatheterized)  hind  leg.  In 
py piglets  2  and  3,  a  second  probe  ( probe  PL)  was  placed  on  the 
M'Catheterized  (left)  leg.  The  protocol  consisted  of  varying  the 
femoral  Svo2  over  the  approximate  range  of  20-95%  by  mod¬ 
ulating  the  volume  fraction  of  oxygen  inspired  by  the  piglet 
f Fign  over  the  range  of  10-100%.  The  oxygenation  cycles 
'pSFmrmed  on  the  three  piglets  are  illustrated  in  Fig.  3.  Each 
cycle  consisted  of  varying  the  Flo,  approximately  every  4-6 
min  through  the  values  of  ~40,  t5,  10,  and  100%  (piglets  1 
and  2 )  or  ~40,  17.5,  15, 12.5,  10,  and  100%  (piglet  3).  We 

performed  two(FTj\  cycles  on  piglet  1,  four  on  piglet  2,  and 
three  on  piglet  3\  For  each  specific  value  of  Fio2,  we  acquired 
about  3,000  optical  data  points  [4  min  X  (60  s/min)/(80 
ms/data  point)]  or  more.  During  cycles  C  and  D  on  piglet  2, 
the  optical  probe  PR  was  slightly  moved  with  respect  to  the 
location  examined  during  cycles  A  and  B,  to  collect  data  on 
two  different  muscle  volumes  during  the  two  cycle  pairs  A-B 
and  C-D.  Optical  probe  PR  always  collected  data  on  the  right 
hind  leg,  whereas  probe  PL  was  placed  on  the  left  hind  leg 
during  cycles  A  and  B  of  piglet  2  and  cycles  A  and  B  of  piglet 
3  (we  did  not  collect  data  with  the  optical  probe  PL  on  piglet 
1,  during  cycles  C-D  on  piglet  2,  and  during  cycle  C  on  piglet 
3).  In  all  three  piglets,  the  invasive  measurement  of  Svo2 
from  a  femoral  vein  blood  sample  was  performed  at  the  end  of 
each  Fio2  interval,  as  shown  in  Fig.  3.  Motion  artifacts  were  ■ 
minimized  in  the  optical  data  by  securing  the  piglet’s  legs  to . 
the  operating  table.  The  protocol  was  approved  by  the  Insti¬ 
tutional  Review  Board  of  the  Massachusetts  General  Hospi¬ 
tal,  where  the  piglet  experiments  were  performed. 

Measurements  on  human  subjects.  We  performed  measure¬ 
ments  on  eight  healthy  human  subjects  (6  men  and  2  women; 
mean  age  of  24.5  yr,  age  range  of  20-35  yr).  The  subjects  sat 
on  a  comfortable  chair  and  rested  for  10-15  min  before  the 
experimental  protocol  was  started.  A  pneumatic  cuff  was 
placed  around  the  right  thigh  of  the  subject  to  later  induce  a 
venous  occlusion  by  inflating  the  cuff  to  a  pressure  of  70 
mmHg.  A  pulse  oximeter  probe  (Nellcor,  N-200)  was  placed 
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AQ:  XI  Fig.  2.  Experimental  arrangement  for  the  piglet  study.  A  breathing  mask  applied  to  the  piglet’s  snout  provided  the 
3-4%  isoflurane  anesthetic  and  was  connected  to  the  oxygen  line  for  variations  in  the  fraction  of  inspired  oxygen 

O(Fl0,).  A  strain-gauge  belt  and  a  pulse  oximeter  monitored  the  respiratory  excursion  and  the  heart  rate, 
0  respectively,  and  their  analog  outputs  were  directed  to  the  auxiliary  inputs  of  the  frequency-domain  tissue 
spectrometer  (ISS,  Champaign,  IL,  model  96208).  One  or  two  optical  probes  (PR  on  the  right  hind  leg  and  PL  on 
the  left  hind  leg)  of  the  tissue  spectrometer  were  used  to  measure  the  near-infrared  tissue  absorption  with  a  time 
resolution  of  80  ms.  The  absorption  oscillations  at  the  respiratory  frequency  were  processed  to  provide  measure¬ 
ment  of  the  venous  02  saturation  (Svfe'-NTRSreBp)  (NIRS  is  near-infrared  spectroscopy).  Invasive  measurements  of 
, — x  the  venous  O2  saturation  (designated  Svo,,-bIaad)  were  obtained  by  gas  analysis  of  venous  blood  samples  collected 
AQ:  8^by  a  femoral  vein  catheter,  aux,  Auxiliary.  i/ot  opt^tz=>'-  'upcrr^-po  ipjx 
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Fig.  3.  Schematic  representation  of 
FfSj  cycles  for  piglet  1  (A),  piglet  2  ( B ), 
and  piglet  3  (C).  •,  Time  at  which  ve- 
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on  the  index  finger  of  the  left  hand.  A  strain-gauge  belt 
(Sleepmate/Newlife  Technologies,  Resp-EZ)  was  placed 
around  the  subject’s  upper  abdomen  to  monitor  the  respira¬ 
tory  excursion.  As  in  the  piglet  experiment,  we  used  the 
analog  outputs  of  the  pulse  oximeter  and  strain  gauge  for 
continuous  coregistration  of  the  physiological  and  near-infra- 
red  data.  Two  optical  probes  were  placed  on  the  right  thigh; 
the  first  probe  (probe  HVM)  was  positioned  on  top  of  a  visible 
superficial  vein  of  the  vastus  medialis  muscle,  and  the  second 
probe  ( probe  HVL)  was  placed  on  the  vastus  lateralis  muscle, 
far  from  visible  superficial  veins.  During  the  measurements, 
we  asked  the  subject  to  breathe  regularly,  following  a  met¬ 
ronome  whose  frequency  was  set  to  the  average  breathing 
rate  of  the  subject  at  rest  (typically  14-15  breaths/min). 
During  the  whole  experiment,  the  subject  was  asked  to 
breathe  at  the  same  frequency  as  the  metronome  pace.  No 
subjects  experienced  any  discomfort  or  difficulties  with  this 
procedure.  The  measurement  protocol  consisted  of  2  min  of 
baseline  (we  acquired  1,535  optical  data  points  at  80  ms/ 
point),  followed  by  40  s  of  venous  occlusion,  and  a  final 
recovery  period  of  a  few  minutes.  A  few  subjects  performed 
an  additional  exercise  routine  to  test  the  effect  of  exercise  on 
the  measured  value  of  Svo2-NIRSresp  on  the  muscle.  The 
exercise  consisted  of  raising  the  right  foot,  voluntarily  con¬ 
tracting  the  leg  muscles  (isometric  contraction),  until  the 
subject  felt  tired.  The  human  study  was  approved  by  the 
Institutional  Review  Board  of  Tufts  University,  where  the 
human  experiments  were  performed;  all  subjects  gave  their 
written,  informed  consent. 

Near-infrared,  data  processing  for  the  measurement  of  SvQ. 
We  used  a  modified  Beer-Lambert  law  approach  (14)  to 
translate  the  temporal  intensity  ratio  collected  at  each  wave¬ 
length  [/(X,  t)//(Xb@),  where  I  is  intensity,  X  is  wavelength, 
and-iisJimej  at  a  distance  of  1.0  cm  from  the  illumination 
point,  into  a  time  variation  in  the  tissue  absorption  [Apa(X,  t). 


where  pa  is  the  tissue  absorption  coefficient].  This  approach 
was  implemented  by  applying  the  following  equation  (14) 


Ap0(X,  f)  =  — In 


/(X,  0) 


U(M)J 


(. 1 ) 


where  Lea  is  the  effective  optical  pathlength  from  the  illumi¬ 
nating  point  to  the  light  collection  point.  We  measured  Lefrby 
quantifying  pa  and  the  reduced  scattering  coefficient  (pi) 
using  the  frequency-domain  multidistance  method  (17).  The 
diffusion-theory  relationship  that  gives  Les  in  terms  of  p*,, 
pj,  and  the  source-detector  separation  (r)  in  a  semi-infinite 
turbid  medium  (where  the  illumination  and  collection  points 
are  at  the  boundary  of  the  turbid  medium)  is  the  following 
(17) 


Lett  ~ 


3p>2 
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More  details  on  this  hybrid  frequency-domain  [to  measure 
Let/X)]  and  continuous  wave  (modified  Beer-Lambert  law) 
approach  are  given  in  Refs.  14,  17,  21,  and  22.  Equation  2 
shows  that  for  typical  values  of  the  near-infrared  pa  and  pi, 
say  p0  =  0.1  cm-1  and  pi  =  10  cm-1,  the  value  of  Leffis  ~5.5 
cm  for  r  =  1  cm.  The  multi-distance  scheme  was  implemented 
by  considering  the  data  collected  by  the  two  fiber  bundles 
located  at  two  different  distances  (1.0  and  2.0  cm)  from  the 
source  fiber  bundle.  At  these  source  detector  distances,  the 
diffusion  regime  of  light  propagation  in  tissues  is  already 
established  (18).  As  an  alternative  to  the  diffusion  equation 
model  to  describe  the  spatial  dependence  of  the  optical  signal, 
empirical  approaches  have  been  proposed  (6).  The  different 
sensitivity  of  the  two  detector  channels  was  accounted  for  by 
a  preliminary  calibration  measurement  on  a  synthetic  tissue¬ 
like  sample.  The  applicability  of  the  initial  calibration  to  the 
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whole  data  set  was  verified  at  the  end  of  each  measurement 
session  by  repositioning  the  optical  probe  on  the  calibration 
sample.  We  typically  reproduced  the  calibration  values  of  the 
block  optical  coefficients  to  within  10%.  In  the  piglet  experi¬ 
ments,  we  updated  the  measured  values  of  Leg  at  each 
wavelength  every  time  the  Fio2  was  changed.  Specifically, 
Leff  was  computed,  according  to  Eq.  2,  from  average  measure¬ 
ments  of  tv  and  pj  over  the  last  80  s  of  each  period  corre¬ 
sponding  to  a  specific  Fio2  value.  For  the  measurements  on 
human  subjects,  we  computed  an  initial  value  of  Les  (at  each 
wavelength)  over  the  first  2  min  of  baseline,  and  we  used  this 
value  for  the  analysis  of  the  data  over  the  whole  measure¬ 
ment  session.  The  long  integration  time  for  the  mean  path- 
length  measurements  (80  s  in  the  piglet  experiment,  120  s  in 
the  human  subjects  measurements)  realized  a  low-pass  filter 
that  minimized  the  time-varying  contributions  from  the  Hb 
oscillations  caused  by  the  arterial  pulsation  and  breathing. 
Furthermore,  the  two-distance  measurement  scheme  for  the 
mean  pathlength  measurement  also  provided  some  level  of 
spatial  averaging.  In  contrast,  the  optical  data  for  the  mea¬ 
surement  of  Svo2  were  acquired  with  an  80-ms  temporal 
resolution  and  with  the  use  of  a  single  source-detector  dis¬ 
tance  (1  cm).  - 

To  measure  the  Svo2,  we  followed  a  two-step  procedure. 
First,  we  computed  the  amplitude  of  the  absorption  oscilla¬ 
tions  at  the  respiratory  frequency  at  each  of  the  eight  wave¬ 
lengths  considered.  Second,  we  fit  the  spectrum  of  the  ex¬ 
perimental  absorption  amplitude  with  the  hemoglobin  ab¬ 
sorption  spectrum.  We  have  used  two  alternative  methods  to 
quantify  the  absorption  oscillations  at  the  respiratory  fre¬ 
quency.  The  first  method  is  based  on  the  fast  Fourier  trans¬ 
form  (FFT)  of  A  \La(t).  The  sum  of  the  amplitudes  of  the  FFT  of 
Ajv  over  the  respiratory  frequency  band  yields  a  measure  of 
the  amplitude  of  the  respiration-induced  absorption  oscilla¬ 
tions.  This  method  assumes  that  the  Fourier  spectrum  of  Ap„ 
clearly  shows  a  discemable  peak  at  the  respiratory  fre¬ 
quency.  The  second  method  is  based  on  a  band-pass  (BP) 
filter  of  Ap,0(t)  and  on  a  modeling  algorithm  (MA)  (sine-wave 
fit).  The  BP  filter  serves  the  purpose  of  isolating  the  absorp¬ 
tion  oscillations  at  the  respiratory  frequency  by  suppressing 
higher  and  lower  frequency  components  in  Apa(£).  The  MA 
consists  of  fitting  a  sine  wave  to  ApQ(BP)  over  each  respira¬ 
tory  cycle.  The  amplitude  of  the  fitted  sine  wave  gives  an 
estimate  of  the  absorption  oscillation  amplitude  at  the  respi- 1 
ratory  frequency.  As  a  result,  the  second  method  (BP  +  MA)/ 
achieves  a  reading  of  Svo2  from  each  individual  respiration 
cycle,  whereas  the  first  method  (FFT)  requires  multiple,  res-! 
piration  cycles  to  produce  a  Svo2  reading.  Both  methods! 
provide  phase  readings  that  can  be  used  to  verify  that  the 
respiration-induced  absorption  oscillations  at  different  wave¬ 
lengths  are  in  phase  with  each  other.  We  indicate  the  Svo2 
measurement  according  to  the  FFT  and  BP  +  MA  methods 
with  Svo2-NIRSraBp(FFT)  and  fSvo2-NJRSroSp(RP),  respec,-/ 

tively.  ^ - ~  " 

In  the  piglet  experiments,  we  evaluated  the  FFT  of  Apa 
over  256  data  points,  corresponding  to  a  time  trace  of  20.5  s, 
to  achieve  reliable  spectra  from  a  number  of  breathing  peri¬ 
ods  (typically  13-16).  Furthermore,  we  averaged  about  800 
successive  FFTs  (each  computed  from  a  data  set  shifted  by 
one  data  point  with  respect  to  the  previous  one),  so  that  the 
total  number  of  data  points  resulting  in  a  single  Svo2  reading 
was  on  the  order  of  1,000,  corresponding  to  a  train  of  data 
80  s  long.  This  80-s-long  data  set  was  chosen  to  be  at  the  end 
of  each  Fio2  period,  and  it  coincides  with  the  80-s  period  over 
which  we  measured  Leg.  In  the  human  subject  experiment, 
we  used  512  points  for  the  FFT  because  the  breathing  fre¬ 
quency  was  lower  (0.22—0.26  Hz)  them  that  of  the  piglets 


(0.6-0.9  Hz)  and  we  wanted  to  have  a  similar  number  of 
breathing  periods.  As  in  the  piglets  experiment,  we  averaged 
the  results  from  multiple  (500-1,000)  successive  FFTs. 

The  spectrum  of  the  amplitude  of  the  absorption  oscilla¬ 
tions  at  the  respiratory  frequency  [Ap£esp(A;)]  was  fitted  with 
a  linear  combination  of  the  Hb02  and  Hb  extinction  spectra, . 
|^j%)A[Hb02]^p  +  2Hb(A,)A[Hb]-^  where/fei^a,) 
are  the  extinction  coefficients  of  HbC>2  EttulHb, 
respectively  (43,  53).  The  fitting  parameters  were  the  ampIF" 
tudes  of  the  oscillatory^-EHbQa3^tA{HfeOd£S^-;And_EHbl 
(A[IJb3r®^p)  at  the  respiratory  frequency.  The  mimrrfkration  of 
the  sum  of  the  squares  of  the  residuals,  i.e.,  2;[A|i^(X;)  - 
Ap.£eep(A;)]2,  yields  a  linear  system  whose  solution  gives\the 
following  best  fit  concentrations  of  amplitude  of  the  oscilla¬ 
tory  [HbC>2]  and  [Hb]  (11)  \ 


to  s 


A[Hb02]resp 


A[Hb]reap  = 


[2i  Ap™p(A;)£Hbo20'-;)]R  EHbQO]  \ 

-  [2;  Ap.„liaP(A;)Em(A,)][S;  £Hb02(^i)EHb(^t)]  \ 

[2;  EHb02(A;)][2;  £[ib(A,)]  -  [2;  £Hb02(A,-)£Hb(Ai)]2  \ 

(3)  6x.c 

[2,  Ap.rp(Ai)£Hb(A;)][2,'  eLo2(A)] 

~  [2,  A^aP(A,)£mo2(\,)][2,  EHb02(ki)£Ht(A;)3 
[2/  eHbo2(>0][2;  EHb(A,)]  -  [2,  EHbo,(Ai)£Hb(A,)]2 


The  oxygen  saturation  of  the  hemoglobin  compartment  oscil¬ 
lating  synchronously  with  respiration  (Svo2-NIRSresP)  is  then 
given  by 


-  A[Hb02r8p 

Sv$  -N.IRSro8p  -  ApIb0jMp  +  A[Hb]r.»P  0 

[2;  Ap-‘p(Xi)£Hb02(Ai)][2i  EHb(A;)] 

-  [2;  Ap.™”p(\i)eHb(\i)][2;  Enb02(A)EHb(\)] 


t* _  [2;  Ap.™8p(\;)eHb(\i)][2;  Enb02(A)EHb(\)] 

[2;  Ap.„esp(\i)EHbo2(ki)]}2i  £Hb(A,)[£Hb(Xi)  -  £Hbo2(A,)]6jic  y 
I  ~  [2i  A|i.™"P(\,)eHb(X;)]  {2,  £Hb02(^i)[£Hb(Xi)  —  EHb02(Ai)]}  * 

/It  is  important  to  note  that  for  the  determination  of  Svo2- 
I  NIRSresP  one  only  needs  to  know  Z,effto  within  a  wavelength- 
I  independent  factor.  In  fact,  Eq.  5  shows  that  a  common, 
j  wavelength-independent  multiplicative  factor  in  Ap.a(A ;)  can- 
'  cels  out  in  the  expression  for  Svo2-NIRSresP.  In  contrast,  the 
wavelength  dependence  of  Leg  is  important  for  the  measure¬ 
ment  of  Svo2  with  our  method,  and  this  is  why  we  have  opted 
to  measure  Leg  at  each  wavelength  using  the  multidistance, 
frequency-domain  technique.  It  is  also  important  to  observe 
\  that  our  method  requires  that  1)  oscillations  of  p*,  at  the 
\  respiratory  frequency  to  be  reliably  attributed  to  hemoglobin 
j  (and  not,  for  instance,  to  motion  artifacts),  2)  the  hemoglobin 
/  concentration  fluctuations  to  result  from  the  volume  oscilla- 
'  tion  of  a  hemoglobin  compartment  rather  than  from  periodic 
fluctuations  in  the  blood  flow,  and  3)  the  fluctuating  hemo¬ 
globin  compartment  responsible  for  the  measured  Ap„  to  be 
mainly  the  venous  compartment.  In  our  measurements,  we 
have  considered  each  one  of  the  three  above  points.  The 
assignment  of  the  absorption  oscillations  to  hemoglobin 
( point  1 )  was  done  by  requiring  that  the  hemoglobin  spec¬ 
trum  fits  the  absorption  data  relatively  well.  To  this  aim,  we 
requested  that  the  average  absolute  value  of  the  relative 
residuals,  defined  as  (2^=  1/1^1  jfApJj^Ai)  -  Apf7°p(A;)l/ 
Ap®*(Xi),  where  N  is  tEenunibeFffwavelimgths  considered,’ 
"be  at  "most  twice  the  experimental  percent  error  in  Ap™Bf>.  We 
also  used  the  standard  deviation  of  the  Svo2-NIRSresp(FFT) 
values  obtained  with  the  800  (piglet  experiment)  or  500- 
1,000  (human  experiment)  successive  FFTs  to  estimate  the 
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error  in  S^2-NIRSresp(FFT).  We  discarded  the  cases  having 
a  standard  deviation  error  in  Svo2-NIRSreap  greater  than 
15%.  The  assignment  of  the  absorption  oscillations  to  volume 
rather  than  blood  flow  fluctuations  (point  2)  is  achieved  by 
verifying  that  the  absorption  oscillations  at  the  eight  wave¬ 
lengths  are  in  phase.  In  fact,  blood  flow  fluctuations  induce 
out-of-phase  oscillations  in  the  [Hb02]  and  [Hb]  (because  of 
the  increased  rates  of  inflow  of  Hb02  and  washout  of  Hb),  as 
opposed  to  the  in-phase  oscillations  of  Hb02  and  Hb  that 
result  from  volume  pulsations.  The  third  point,  namely  the 
requirement  that  the  absorption  oscillations  at  the  respira¬ 
tory  frequency  are  representative  of  venous  blood,  is  investi¬ 
gated  by  1)  comparing  the  Svo2-NIRS  from  the  respiratory 
hemoglobin  oscillations  (Svo2-NIRSresP)  with  the  correspond¬ 
ing  values  measured  by  gas  analysis  of  Svo2-blood  (piglet 
experiments)  or  by  the  NIRS  venous  occlusion  method  (Svo2- 
NIRSvo)  (human  subject  experiments),  2)  studying  the  effect 
on  the  [Hb]  and  [Hb02]  oscillations  at  the  respiratory  fre¬ 
quency  of  a  venous  occlusion  induced  between  the  lungs  and 
the  peripheral  measurement  area  (the  thigh  muscles  in  hu¬ 
man  subject  experiments),  and  3)  by  recording  the  effect  of 
muscle  exercise  on  the  near-infrared  measurements  of  Svo2 
[Svo2-NIRSresP(BP)]  in  human  subjects.  f  “!  *  ? 

RESULTS 

Figure  4  reports  average  spectra  of  Leff  measured  for 
a  source-detector  separation  of  1  cm.  Figure  4A  refers 
to  piglet  measurements  conducted  at  two  different  val¬ 
ues  of  F$2,  whereas  Fig.  4 B  refers  to  human  measure¬ 
ments  with  probes  HVM  and  HVL.  The  error  bars  in 


Fig.  4.  Near-infrared  spectra  of  the  effective  optical  pathlength  (Le n) 
measured  on  the  piglet’s  leg  (A)  and  a  human  thigh  muscle  (B)  for  a 
source-detector  separation  (r)  of  1  cm.  In  A,  different  symbols  refer  to 
2  different  values  of  FlQ2.  In  B,  different  symbols  refer  to  2  different 
thigh  muscles  (vastus  medialis  for  probe  HVM  and  vastus  lateralis 
for  probe  HVL).  The  lines  join  the  points  as  an  aid  to  the  eye. 


Fig.  4  represent  the  standard  deviations  over  multiple 
measurements  (multiple  Fio2  cycles  and  piglets  for  Fig. 

4 A  and  multiple  subjects  for  Fig.  4B). 

In  the  piglet  experiment,  we  discarded  11  (from  a 
total  of  67)  Sv@2-NIRSreSp(FFT)  measurements  because 
the  standard  deviation  over  800  FFTs  exceeded  15%. 

These  discarded  Svo2-NIRSresp(FFT)  readings  occurred 
as  follows:  one  (of  8)  in  piglet  1,  two  (of  26)  in  piglet  2, 
and  eight  (of  33)  in  piglet  3.  One  discarded  reading  was 
assigned  to  motion  artifacts,  whereas  the  other  ten 
discarded  measurements  all  occurred  at  low-Fin.  val-  . 
ues  (10-17.5%)  corresponding  tc(  Svcybloocpralues  of  &1-OC 
20-50%.  We  were  not  able  to  apply  the  Hpmethod  to 
piglet  2  and  to  the  Fio2  cycles  A  and  B  of  piglet  3 
because  of  irregular  absorption  oscillation  waveforms 
that  were  not  reliably  processed  by  the  HP  +  MA#aQ:4J 

approach.  ' . 

•••  Figure  5  shows  typical  temporal  traces  of  the  relative  fs‘3*-« 
[Hb02]  and  [Hb]  measured  on  the  piglet’s  leg  (with 
optical  probe  PR)  (Fig.  5 A)  and  on  the  human  vastus  ° ' 6 
medialis  muscle  at  rest  (Fig.  5 B)  and  during  venous 
occlusion  on  the  upper  thigh  (optical  probe  HVM)  (Fig. 

5 C).  The  temporal  traces  of  [Hb]  and  [Hb02]  are  ob¬ 
tained  by  fitting  the  measured  spectrum  of  Ap.a(\,f) 

(whose  value  at  each  wavelength  was  obtained  from 
Eq.  2)  with  a  linear  combination  of  the  Hb02  and  Hb 
extinction  spectra.  This  procedure  results  in  the  appli¬ 
cation  of  Eqs.  3  and  4  without  the  superscript  “resp”  on 
A(jia,  A[Hb02],  and  A[Hb],  Two  oscillatory  components 
are  clearly  visible  in  the  relative  [Hb02]  and  [Hb] 
traces  of  Fig.  5A:  the  first  one,  associated  with  the 
heartbeat  (as  shown  by  the  pulse  oximeter  data;  top 
trace  in  Fig.  5)  is  at  a  frequency  of  —2.5  Hz,  whereas 
the  second  one,  associated  with  respiration  (as  shown 
by  the  strain  gauge  signal;  second  trace  from  the  top  in 
Fig.  5),  is  at  a  frequency  of  —0.65  Hz.  Only  the  latter 
oscillatory  component  (at  a  frequency  of  —0.23  Hz  in 
human  subjects)  is  clearly  visible  in  Fig.  5B,  whereas 
neither  is  present  in  Fig.  5C.  Figure  5,  B  and  C,  shows 
additional  low-frequency  oscillations  associated  with 
changes  in  blood  pressure  and  heart  rate.  We  observe 
that  the  strain-gauge  signal  (second  trace  from  the  top 
in  Fig.  5)  increases  during  inspiration  and  decreases 
during  expiration.  The  BP  filter  described  in  the  pre¬ 
vious  section  aims  at  isolating  the  oscillatory  compo¬ 
nent  at  the  respiratory  frequency  by  filtering  out 
higher  and  lower  frequency  components.  The  relative 
[Hb02]  and  [Hb]  traces  after  BP  filtering  are  shown  in 
Fig.  5,  bottom.  In  the  case  reported  in  Fig.  5,  which  is 
representative  of  the  results  reported  in  this  article  for 
Svo2-NIRSresp,  the  oscillatory  components  of  [Hb02] 
and  [Hb]  at  the  respiratory  frequency  are  in  phase  with 
each  other  and  disappear  during  venous  occlusion. 

Figure  6  illustrates  representative  Aix'eBp  spectra  F6 
measured  on  the  piglet’s  leg  ( probe  PR)  (Fig.  6A),a«4  on  - 
the  human  vastus  medialis  muscle  at  rest  (Fig.  6B),and  (y — 
during  venous  occlusion  on  the  upper  thigh  (probe 
HVM)  (Fig.  6C).  They-axis  of  each  panel  of  Fig.  6  refers 
to  the  values  of  A|xjjesp  obtained  with  the  BP  filter 
method.  The  values  of  A|i£esp  computed  with  the  FFT 
method  are  normalized  by  a  wavelength-independent 


J  Appl  Physiol  •  VOL  91  •  DECEMBER  2001  •  www.jap.org 


tapraid4/dg-japp/dg-japp/dg1 201/dgl  1 69d01  o  thoenyj  S=9  9/13/01  18:36  MS:  A0506-1  |  Ini:  12e/rgh  [  lnput-DCT-123 
MEASURING  VENOUS  SATURATION  WITH  NEAR-INFRARED  SPIROXIMETRY 


piglet  (FiO2=40%) 


B 


human  subject  at  rest 


human  subject 
during  venous  occlusion 


heartbeat 

breathing 

signal 

(a.u.) 

1 

d 

.2 

1 

T3 

I 

£ 

§ 

a 

Ch 

o 

o 

.2 

o 

w 

i 

< 

JS 

<u 

s 

4. 

.5 

03 

1 

T 
/— s 

Oh 

CO 

\AAAAAA  /WWW  vVVVWV 


10  15  20 

time  (sec) 


Fig.  5.  Representative  traces  of  the  relative  HbC>2  concentrations  ([HbOzD  and  Hb  concentrations  011b])  measured 
on  the  piglet’s  leg  (with  optical  probe  PR)  (A)  and  on  the  human  vastus  mcdialis  muscle  at  rest  (B)  and  during 
venous  occlusion  on  the  upper  thigh  (optical  probe  HVM)  (C).  Bottom  panels  report  the  [Hb]  traces  after  processing 
with  the  digital  band-pass  filter  [band  pass  (BP)  +  modeling  algorithm  (MA)[  designed  to  isolate  the  oscillations 
at  the  respiratory  frequency.  Top  trace  represents  the  piglet’s  heartbeat  monitored  by  the  pulse  oximeter.  Second 
trace  from  the  top  is  the  strain  gauge  signal  that  monitors  the  respiratory  excursion.  The  strain  gauge  signal 
increases  during  inspiration  and  decreases  during  expiration,  a.u.,  Arbitrary  units. 


(AQ:  5)  factor  to  match  the  BP  value  of  Ap,™sp  at  636  nm.  The 
relatively  high  value  of  (e^M  during  venous  occlusion 
(Fig.  6C)  is  an  indication  oFthe  poor  fit,  which  in  turn 
results  from  the  lack  of  hemoglobin  oscillations  at  the 
5  respiratory  frequency  (see  Fig.  5 C,  bottom,  and  the 

©y-axis  values  of  Fig.  5 C  compared  with  those  of  Fig. 
5 B).  Figure  6  also  shows  the  best  fit  of  the  hemoglobin 
absorption  spectrum  to  the  BP  Ap£esp  and  to  the  FFT 
Ap£esp.  The  best-fit  hemoglobin  spectra  represent  the 
oxygen  saturation  of  hemoglobin,  as  illustrated  in  Fig. 
1.  The  value  of  Svo2-NIRSreSp  is  given  by  Eq.  5. 

F7  Figure  7  compares  the  measurements  of  Svo  - 
NIRSresp(BP),  Svo2-NIRSresp(FFT),  and  Svo2-blood  dur¬ 
ing  cycle  A  of  piglet  1  and  during  cycle  C  of  piglet  3.  The 
Svo2-NIRSreSp(BP)  traces  reported  in  Fig.  7  were  ob¬ 
tained  by  performing  a  running  average  of  the  breath- 
to-breath  values  obtained  with  the  BP  method.  In  Fig. 
7,  the  averaging  procedure  consists  of  a  5-point  (in  Fig. 
7A)  or  15-point  (in  Fig.  IB)  running  average.  The 
assessment  of  the  agreement  between  the  measure¬ 
ments  of  Svo2-NIRSresp(FFT)  andJ3vg^BloocHn  the  full 
piglet  study  is  carried  out  accofdihg4-crtheprocedure 
described  by  Bland  and  Altman  (5).  Figure  8A  plots  the 


results  of  the  NIRS  method  based  on  the  respiratory 
oscillations  of  the  tissue  absorption  against  the  inva¬ 
sive  measurement  of  Svo2-blood.  The  shape  of  the  sym¬ 
bols  in  Fig.  8A  indicates  the  piglet  number,  whereas 
the  type  of  fill  indicates  the  location  of  the  NIRS  mea¬ 
surement.  The  range  of  Svo2-blood  values  considered  in 
this  study  is  —20-95%.  The  error  bars  in  Fig.  8A  are 
the  standard  deviations  computed  from  the  results  of 
—800  successive  FFTs  (as  described  in  materials  and 
methods).  Figure  8 B  displays  the  difference  between 
the  two  readings  vs.  their  average,  and  it  quantifies  the 
discrepancy  between  the  two  methods  and  the  possible 
dependence  of  such  a  difference  on  the  level  of  Svo2. 
The  mean-difference  between  Svo2-NIRSresp  and  Svo2- 
blood  over  the  fufloxygenation  range  considered  in  this 
study  is  1.0%  (a  measurement  of  the  bias  of  the  Svo2- 
NIRSresp  measurement),  and  the  standard  deviation  pf 
the  difference  is  5.8%.  Figure  8 B  does  not  show  any 
striking  dependence  of  the  difference  on  the  mean.  If 
we  take  the  values  of  mean  difference  ±  2  SD  as  the 
limits  of  agreement  of  the  two  methods  (5),  we  get  an 
estimate  of  the  maximal  discrepancies  between  Svo2- 
NIRSreSp(FFT)  and  Svo2-blood  of  -10.6%  and  +12.6%. 


SD 
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piglet  (FiO2=40%) 
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human  subject  at  rest 


human  subject  during  venous  occlusion 


Fig.  6.  Representative  change  in  respiratory  tissue  absorption  coef¬ 
ficient  (4|i'“p)  spectra  measured  with  the  BP  and  fast  Fourier 
transform  (FFT)  methods  on  the  piglet’s  leg  (probe  PR)  (A)  and  on  the 
human  vastus  medialis  muscle  at  rest  (B)  and  during  venous  occlu¬ 
sion  on  the  upper  thigh  (probe  HVM)  (C).  The  experimental  8-point 
A|a£“p  spectra  were  fitted  with  the  hemoglobin  absorption  spectrum 
(with  the  hemoglobin  concentration  a*Aee*ewtite<las  fitting  param¬ 
eters).  The  values  of  efit  (defined  in  the  text)  for  the  BP  and  FFT 

i  Spectra  give  a  measure  of  the  quality  of  the  fit. 

- - 

j  ,J>V  In  the  human  experiment,  we  found  that  the  NIRS 
values  of  Svo2  measured  with  probe  HVM  (placed  on 
top  of  a  visible  vein)  were  typically  smaller  than  those 
measured  with  probe  HVL  (placed  far  from  any  visible 
vein).  Furthermore,  the  amplitude  of  the  oscillatory 
absorption  at  the  respiration  (heartbeat)  frequency 
was  typically  greater  (smaller)  for  the  data  collected 
with  probe  HVM  than  with  probe  HVL.  Of  the  16 
Svo2-NIRSresp  measurements  (8  subjects,  2  locations), 


we  discarded  only  2  measurements  (because  of  a  value 
of  2flt  greater  than  twice  the  error  in  A|x£esp),  both 
collected  with  probe  HVL.  Figure  9  compares  the  Svo2-  F9 
NIRSresp  values  measured  in  the  human  subjects  at 
rest  in  which  the  FFT  method  and  the  BP  filtering 
approach  were  used.  Figure  9A  shows  the  good  agree¬ 
ment  of  the  two  measurements,  and  Fig.  9 B  quantifies 
the  average  difference  (0.9%)  and  the  maximum  dis¬ 
crepancies  of  -5.1  and  +6.9%,  as  given  by  the  mean  ± 
twice  the  standard  deviation  of  the  differences.  Figure 
10  reports  a  similar  comparison  between  Svo2-NTRSre-  Fio 
sp(FFT)  and  Svo2-NIRSV0.  As  described  by  Yoxall  and 
Weindling  (56),  under  the  assumption  that  a  venous 
occlusion  induces  an  initial  increase  in  the  venous 
blood  volume,  Svo2-NIRSvo  is  given  by  [Hb02]o/[Hb 
-T]o,  where  the  dots  indicate  a  time  derivative  and  the 
subscript  0  indicates  the  initial  time  that  immediately 
follows  the  onset  of  venous  occlusion.  The  agreement 
between  Svo2-NIRSreBp(FFT)  and  Svo2-NIRSvo  is  good, 
with  an  average  deviation  of  0.8%  and  maximum  dis¬ 
crepancies  of  -4.2  and  +5.8%.  Two  horizontal  lines  in 
Figs.  85  and  95  indicate  the  range  given  by  the  mean 
difference  ±  2  SD.  The  maximum  discrepancy  among 
Svo2-NIRSresp(FFT),  Svo2-NIRSresP(BP),  and  Svo2- 
NIRSV0  is  less  than  the  maximum  deviation  between 
Svo2-NIRSreBp(FFT)  and  Svo2-blood  found  in  piglets 
(see  Fig.  86). 

The  effect  of  muscle  exercise  on  the  measurement  of 
Svo2-NIRSr0sP(BP)  on  top  of  a  visible  superficial  vein 
(probe  HVM)  is  illustrated  in  Fig.  11.  Although  Sao2  Fii 
(measured  with  a  pulse  oximeter)  is  unaffected  by  the 
exercise,  Svo2-NrRSre8p(BP)  shows  a  significant  postex¬ 
ercise  decrease  from  a  baseline  value  of  75-78%  down 
to  a  minimum  value  of  —54%.  The  recovery  to  the 
baseline  value  of  Svo2-NIRSreap(BP)  occurs  after  —30  s. 

By  using  the  BP  approach,  we  could  monitor  Svo2- 
NIRSreBP  at  every  breathing  period,  i.e.,  every  —5  s, 
thus  achieving  a  real-time  monitoring  of  Svo2.  We 
observe  that  we  could  not  obtain  meaningful  measure¬ 
ments  of  Svo2-NIRSresp  during  exercise  because  of  mo¬ 
tion  artifacts. 

DISCUSSION 

Various  methods  for  measuring  Svo2 ■  The  method 
presented  in  this  article  to  measure  Svo2  from  the 
near-infrared  absorption  oscillations  at  the  respiratory 
frequency  (spiroximetry)  can  be  implemented  by  using 
a  FFT  or  a  digital  BP  filter  in  conjunction  with  a  MA. 

We  have  indicated  the  measurements  of  Svo2  obtained 
with  these  two  approaches  with  the  notations  Svo2- 
NIRSresP(FFT)  and  Svo2-NIRSresp(BP),  respectively. 

An  alternative  method  for  measuring  Svo2  with  NIRS 
is  based  on  a  previously  described  venous  occlusion 
protocol  (13,  39,  55,  56).  We  have  identified  the  results 
of  this  measurement  procedure  with  the  notation  Svo2- 
NIRSvo-  In  the  human  study,  the  NIRS  measurements 
were  conducted  at  two  locations  on  the  thigh.  One 
location  was  on  top  of  a  visible  superficial  vein  of  the 
vastus  medialis  muscle  (probe  HVM),  and  the  second 
location  was  far  from  visible  superficial  veins  on  the 
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time  (min)  time  (min) 


Fig.  7.  Comparison  between  the  con¬ 
tinuous  measurement  of  Svo2-NIRS„,»p 
(BP)  and  the  discontinuous  measure¬ 
ments  of  Svo2-NIRSre«p(FFT)  and  Svo2- 
blood.  A  refers  to  cycle  A  of  piglet  1, 
whereas  B  refers  to  cycle  C  of  piglet  3. 
The  values  of  Flo?  (%,  left  y-axes)  dur¬ 
ing  the  experiment  are  indicated  by 
the  shaded  profiles. 


OSv02-BL00D  H  Sv02-NIRS  (FFT)  — Sv02-NIRS  (BP) 


vastus  lateralis  muscle  {probe  HVL).  Finally,  the  inva¬ 
sive  measurement  of  Svo2  performed  by  the  gas  anal¬ 
ysis  of  venous  blood  samples  is  indicated  with  Svo2- 
blood.  In  this  section,  we  discuss  the  different  features 
of  these  measurements  of  Svo2,  and  the  comparison  of 
their  results,  as  reported  in  Figs.  7-9.  S.  ':f 
The  FFT  and  BP  filter  approaches  to  near-infrared 
spiroximetry.  The  major  advantage  of  the  BP  approach 
is  that  it  allows  for  a  real-time  measurement  of  Svo2  by 
providing  a  reading  of  Svo2-NIRSresp(BP)  at  every  res¬ 
piration  cycle.  Consequently,  this  method  is  particu¬ 
larly  effective  during  transients,  as  illustrated  by  the 
recovery  of  the  Svo2-NIRSresp(BPj  traces  corresponding 
to  the  sudden  increase  of  Fio2  to  100%  in  piglets  (see 
Fig.  6,  A  and  B),  or  to  the  end  of  the  exercise  period  in 
human  subjects  (see  Fig.  11).  On  the  other  hand,  the 
BP  filter  +  MA  method  is  susceptible  to  fluctuations  in 
the  respiratory  frequency  and  to  irregular  respiration 
patterns.  This  accounts  for  the  fact  that  we  did  not  get 
reliable  readings  of  Svo2-NIRSreSp(BP)  in  piglet  2  and 
in  Fio2  cycles  A  and  B  of  piglet  3.  The  FFT  method  was 
more  robust,  producing  reliable  readings  in  56  of  67 
cases  (84%)  in  the  piglet  study  and  in  14  of  16  cases 
(87%)  in  the  human  study.  It  is  important  to  observe 
that  10  of  the  11  discarded  readings  in  piglets  occurred 


at  low-Svo2-NIRS-blood  values  (20-50%),  and  one  was 
assigned  to  motion  artifacts.  Both  discarded  readings 
in  the  human  study  were  collected  with  probe  HVL, 
which  was  placed  far  from  visible  veins.  Therefore,  we 
have  found  indications  that  the  measurement  of  Svo2- 
NIRSresp(FFT)  is  particularly  robust  at  Svo2  values 
>50%  (in  piglets)  and  when  the  optical  probe  is  placed 
on  top  of  a  visible  superficial  vein  (in  human  subjects). 
Although  the  FFT  method,  which  is  based  on  a  mea¬ 
surement  of  the  integrated  peak  at  the  respiratory 
frequency,  is  less  sensitive  than  the  BP  method  to 
irregular  respiration  patterns,  it  is  not  applicable  dur¬ 
ing  transients.  In  fact,  we  did  not  obtain  reliable  read¬ 
ings  of  Svo2  when  the  time  frame  used  to  compute 
Svo2-NIRSr08p(FFT)  (80  s  in  piglets,  80-120  s  in  human 
subjects)  included  significant  changes  in  the  Svo2. 
When  both  the  FFT  and  the  BP  methods  can  be  ap¬ 
plied,  they  provide  Svo2-NIRSre3p  measurements  that 
are  in  excellent  agreement,  as  shown  in  Figs.  6  and  9. 
The  differences  between  the  two  measurements  (stan¬ 
dard  deviation  of  3.0%)  are  comparable  with  measure¬ 
ment  errors  and  significantly  less  than  the  maximum 
deviation  between  Svo2-NIRSre8p(FFT)  and  Svo2-blood 
(approximately  ±10%)  observed  in  the  piglet  study  (see 
Fig.  8B). 


SvOrBLOOD  (%)  ([SvOrNIRSres?(FFT)]+[SvOrBLOOD])/2  (%) 


Fig.  8.  Comparison  of  Svo2-NIRSrMp 
(FFT)  and  Svo2-blood  in  the  piglet  study. 
The  shape  of  the  symbols  refer  to  the 
piglet  (circles,  piglet  1 ;  squares,  piglet  2; 
triangles,  piglet  3),  whereas  the  filling 
indicates  the  measurement  side  (filled 
symbols,  right  leg,  i.e.,  probe  PR;  open 
symbols,  left  leg,  i.e.,  probe  PL).  A:  Svo2- 
NIRSre.p(FFT)  is  plotted  vs.  Svo2-blood. 
B:  difference  is  plotted  vs.  the  average  of 
the  2  measurements.  Two  horizontal 
lines  indicate  the  range  given  by  mean 
difference  ±  2  SD  (SD  is  the  standard 
deviation  of  the  difference  between  the  2 
measurements). 


.  PR  piglet#! 


OpL  A  PL 

■PRpiglet#2  A  PRpiglet  #3 
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Fig.  9.  Comparison  of  Svo2-NIRSr«P 
(FFT)  and  Svo2-NTRSr<,sp(BP)  in  the 
human  study.  •,  Vastus  medialis  mus¬ 
cle;  i.e.,  probe  HVM.  O,  Vastus  lateralis 
muscle;  i.e.,  probe  HVL.  Probe  HVM 
was  place  on  top  of  a  visible  superficial 
vein,  whereas  probe  HVL  was  far  from 
visible  veins.  A:  Svo2-NIRSreaP(FFT)  is 
plotted  vs.  Svo2-NTRSreap(BP).  B:  dif¬ 
ference  is  plotted  vs.  the  average  of  the 
2  measurements.  Two  horizontal  lines 
in  B  indicate  the  range  given  by  mean 
difference  ±  2  SD. 

•  <S.~£T  beSp.luaoc4=  _  •  probe  HVM  O  probe  HVL 

^rb-%-  osc  \  o—  tva  ccc  p^»Tt  C 

Mp.asn.re.mp.nt.st  of  Him  rNIRSresP  and  Svo2-NIRSvo.  which  would  indicate  a  blood  flow  rather  than  vol- 
Both  of  these  NIRSfmethoHs  WHIMS  UMJ  the  SVo2  ume  oscillations,  thus  rendering  the  NIRSresp 
(NIRSresp  and  NIRSvo)  rely  on  a  change  in  the  vol-  method  inapplicable.  In  our  human  study,  we  found 
umeTraction  of  venous  blood  in  the  tissue.  The  two  an  excellent  agreement  between  Svo2-NIRSresp(FFT) 
major  differences  between  the  two  methods  are  as  and  Svo2-NIRSVO)  with  a  maximum  deviation  on  the 
follows.  1)  TheTilftS^o  method  requires  an  external  order  of  ±4-5%  (see  Fig.  10). 

perturbation  consisting  of  a , pneumatic-cuff-induced  Optical  probes  PR,  PL,  HVM,  and  HVL.  In  the  piglet 
venous  occlusion,  whereas  the  $J^Sreap  method  is  study,  we  have  found  no  significant  difference  be- 
only  based  on  the  intrinsic  blood  pressure  oscilla-  tween  the  Svo  -NIRSresp(FFT)  data  collected  with 
tions  induced  by  normal  respiration  and  can  be  ap-  probes  PR  (on  the  right  leg)  and  PL  (on  the  left  leg, 
plied  continuously.  2)  The  $2ttSV0  method  can  be  where  the  venous  catheter  was  inserted)  (see  Fig. 
applied  only  to  limbs,  whereas  the  lflR§rcBp  method  8A).  This  result  indicates  that  noninvasive  measure- 
can,  in  principle,  be  applied  to  any  tissue  and  in  ments  of  Svo2  on  one  leg  can  be  meaningfully  corn- 
particular  to  the  brain,  as  already  shown  by  Wolf  et  pared  with  invasive  measurements  of  Svo2  on  the 
al.  (52).  However,  we  stress  that  it  is  always  impor-  other  leg.  In  the  human  study,  we  found  some  differ- 
tant  to  verify  that  the  [Hb02]  and  [Hb]  oscillate  in  ences  between  the  Svo2-NIRS  measurements  with 
phase  at  the  respiratory  frequency  for  the  NIRSresp  probe  HVM  (placed  on  top  of  a  visible  superficial  vein 
method  to  provide  reliable  measurements  of  Svo2.  in  the  vastus  medialis  muscle)  and  with  probe  HVL 
For  instance,  Elwell  et  al.  (16)  reported  out-of-phase  (placed  far  from  visible  veins  on  the  vastus  lateralis 
oscillations  of  [Hb]  and  [Hb02]  in  the  human  brain,  muscle),  As  shown  in  Figs.  8  and  9,  the  Svo„-NIRS 


Fig.  10.  Comparison  of  Svo2-NIRSroBp 
(FFT)  and  Svo2-NIRSvo  (venous  occlu¬ 
sion)  in  the  human  study.  •,  Vastus 
medialis  muscle;  i.e.,  probe  HVM.  o, 
Vastus  lateralis  muscle;  i.e.,  probe 
HVL.  Probe  HVM  was  place  on  top  of  a 
visible  superficial  vein,  whereas  probe 
HVL  was  far  from  visible  veins.  A:  Svo2- 
NIRSrasp(FFT)  is  plotted  vs.  Svo2- 
NIRSvo.  B:  difference  is  plotted  vs.  the 
average  of  the  2  measurements.  Two 
horizontal  lines  in  B  indicate  the  range 
given  by  mean  difference  ±  2  SD. 
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exercise 


Fig.  11.  Continuous  measurement  of  Svo2-NIRSr<u>p(BP)  with  optical 
probe  HVM  (vastus  medialis  muscle,  on  top  of  a  visible  superficial 
vein)  on  a  healthy  human  subject  during  baseline  and  after  isometric 
muscle  exercise  (recovery). 

,<**■  i  f 

readings  (with  both  the  NIRSresp  and  NIRSV0 
method)  of  probe  HVM  (see  Figs.  8  and  9)  were 
typically  smaller  than  the  readings  of  probe  HVL 
(see  Figs.  8  and  9).  We  assign  this  result  to  a  partial 
contribution  from  the  capillary  and/or  arterial  com¬ 
partments  picked  up  by  probe  HVL.  In  fact,  although 
the  optical  data  from  probe  HVM  shown  in  Fig.  4,  B 
and  C,  do  not  show  any  visible  contribution  from  the 
arterial  pulsation,  data  from  probe  HVL  (not  shown) 
do  contain  pulsatile  components  at  the  heartbeat 
frequency.  As  a  result,  we  believe  that  the  optical 
probe  should  be  placed  on  top  of  visible  superficial 
veins  for  more  accurate  readings  of  Svo2-NIRSreap  on 
human  subjects.  We  believe  that  the  reason  that 
Svo2-NIRSrosp  readings  in  the  piglet  study  were  in 
close  agreement  with  the  invasive  measurement  of 
Svo2,  despite  the  evident  arterial  pulsation  in  Fig. 
5A,  is  related  to  the  smaller  extent  of  respiratory 
sinus  arrhythmia  in  piglets  with  respect  to  humans. 
In  fact,  respiratory  sinus  arrhythmia  is  the  main 
origin  of  the  arterial  oscillations  at  the  respiratory 
frequency  (38).  The  larger  role  played  by  respiratory 
sinus  arrhythmia  in  human  subjects  with  respect  to 
piglets  will  probably  require  a  more  careful  interpre¬ 
tation  of  the  optical  data  for  spiroximetry.  However, 
the  results  of  Fig.  11  show  the  practical  applicability 
of  spiroximetry  to  human  subjects,  so  that  we  do  not 
expect  respiratory  sinus  arrhythmia  to  introduce  an 
intrinsic  limitation  of  the  method. 

Noninvasive  vs.  invasive  measurements  of  Svo2-  The 
comparison  between  Svo2-NIRSreap(FFT)  and  Svo2- 
hIood)in  the  piglet  study  shows  a  maximum  deviation 
range  of  -10.6%  to  +12.6%.  The  local  character  of 
the  Svo2  (as  opposed  to  the  systemic  nature  of  the 
Sao2)  requires  some  caution  in  the  comparison  of 
invasive  (Svo2-blood)  and  noninvasive  (Svo2-NIRS) 
measurements  of  Svo2.  In  fact,  in  our  piglet  study, 
Svo2-blood  was  measured  from  blood  samples  drawn 
from  the  femoral  vein,  whereas  Svo2-NIRSresp  was 
measured  with  an  optical  probe  placed  on  the  leg 


muscle.  It  is  likely  that  the  NIRS  oscillatory  signal 
(at  the  respiratory  frequency)  is  not  just  representa¬ 
tive  of  the  femoral  vein  and  may  therefore  be  indic¬ 
ative  of  the  oxygen  consumption  at  different  tissue 
areas  than  those  affecting  the  femoral  vein  satura¬ 
tion.  This  fact  may  not  lead  to  significant  differences 
under  rest  conditions,  but  it  may  be  important  un¬ 
der  stress.  Although  we  found  a  good  agreement 
between  Svo2-NIRSreap(FFT)  and  Svo2-blood  over 
the  whole  range  of  Fio2  values  considered  (see  Fig.  8), 
we  observed  a  meaningfully  greater  standard  devia¬ 
tion  of  the  differences  over  the  20—55%  Svo2-blood 
range  (SD  =  7.8%)  than  in  the  55-95%  range  (SD  = 
3.6%). 

Effect  of  muscle  exercise  on  S v OfHIRSresp(BP).  The 
result  reported  in  Fig.  11  serves  the  purpose  of  further 
illustrating  the  potential  of  the  Svo2-NIRSre8p(BP) 
measurement  approach.  In  fact,  Fig.  11  shows  the 
feasibility  of  monitoring  the  Svo2  in  real  time  on  a 
breath-to:breath  basis  (one  data  point  every  4-5  s). 
Furthermore,  the  baseline  Svo2-NIRSresp(BP)  value  of 
75-78%  and  the  exercise-induced  drop  indicate  the 
venous  origin  of  the  saturation  measurement,  since  the 
Sao2  measurement  provided  by  the  pulse  oximeter 
stayed  constant  at  98  ±  1%  for  the  whole  measurement 
period.  On  the  other  hand,  Kg.  11  reports  only  one 
representative  case,  and  more  studies  are  required  to 
quantify  the  effect  of  muscle  exercise  on  the  measure¬ 
ment  of  Svo2-NIRSre8p.  'j 

In  conclusion,  we  have  presented  a  noninvasive 
approach  to  measure  the  Svo2  in  tissues  from  the 
near-infrared  spectrum  of  the  amplitude  of  respira¬ 
tion-induced  absorption  oscillations.  We  have  imple¬ 
mented  this  approach,  which  we  call  near-infrared 
spiroximetry,  by  processing  the  optical  data  with  a 
FFT  method  or  with  a  digital  BP  filter  in  conjunction 
with  a  MA.  More  sophisticated  data  processing 
schemes  may  further  improve  the  effectiveness  and 
the  range  of  applicability  of  spiroximetry.  The  Svo2 
measurements  reported  in  this  article  complement 
previously  ^demonstrated  NIRS  measurements  of 
Sto2  (8,  27)  and  Sao2  (21).  Therefore,  our  results  may 
lead  to  the  design  of  a  noninvasive  optical  instru¬ 
ment  capable  of  providing  simultaneous  and  real¬ 
time  measurements  of  local  Sao2,  Sto2,  and  Svo2. 
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7.2.1  Lasers  and  arc  lamps 

The  emission  of  continuous-wave  lasers  and  arc  lamps  can  be  modulated  using 
devices  based  on  the  electro-optical  (Pockels  cells)4  or  acousto-optical5,6  effect. 
A  Pockels  cell  is  a  birefringent  crystal  whose  indices  of  refraction  can  be  varied  by 
applying  an  electric  field.  The  application  of  a  time-varying  voltage  to  the  Pockels 
cell  modulates  the  relative  phase  delay  of  the  light  components  polarized  along  the 
two  principal  axes  of  the  cell.  If  this  relative  phase  delay,  or  retardation,  oscillates 
between  0  and  n  when  the  modulated  Pockels  cell  is  sandwiched  between  two 
crossed  linear  polarizers,  each  at  an  angle  of  45  deg  with  respect  to  the  principal 
axes  of  the  cell,  one  achieves  an  intensity  modulator.  In  fact,  no  light  is  transmitted 
when  the  retardation  is  0,  while  all  light  is  transmitted  when  the  retardation  is  7t.  An 
acousto-optic  modulator  is  a  material  that  uses  the  piezoelectric  and  the  photoelas¬ 
tic  effects  to  convert  an  oscillating  electric  field  into  mechanical  vibrations,  which 
in  turn  induce  a  spatially  dependent  index  of  refraction.  When  a  standing  acoustic 
wave  is  established,  the  acousto-optic  crystal  behaves  as  an  oscillating  refractive 
index  grating  that  modulates  the  transmitted  light  by  time-varying  diffraction.  Both 
electro-optic  and  acousto-optic  devices  require  the  light  beam  to  be  collimated.  In 
the  case  of  arc  lamps,  appropriate  collimation  optics  are  required.  Pockels  cells 
provide  effective  modulations  up  to  about  500  MHz,  and  acousto-optic  modulators 
up  to  about  300  MHz. 

The  wavelength  of  the  laser  is  chosen  on  the  basis  of  the  requirements  of 
the  particular  application  (absorption  band  of  a  chromophore,  optimal  penetration 
depth  in  tissues,  etc.).  Examples  of  externally  modulated  cw  lasers  suitable  for  op¬ 
tical  studies  of  tissues  include  the  krypton  ion  (647  nm)  and  He-Ne  (633  nm)  lasers. 
Dye  lasers  pumped  by  either  argon  or  krypton  afford  continuous  tunability  over  a 
wide  spectral  range  that  covers  the  whole  visible  band.  Arc  lamps  (Xe,  Xe-Hg,  etc.) 
provide  continuous  spectral  emission  from  the  UV  (230  nm)  to  the  near-infrared 
(1100  nm).  Therefore,  they  are  ideal  sources  for  spectroscopic  studies  when  a  wide 
and  continuous  spectral  range  is  required. 

7.2.2  Pulsed  sources 

It  is  possible  to  achieve  a  large  modulation  bandwidth  by  exploiting  the  harmonic 
content  of  pulsed  sources  with  high  repetition  rates.  These  sources  can  be  either 
mode-locked  pulsed  lasers  (Nd:  YAG,  Ti: sapphire,  dye  lasers,  etc.)7  or  synchrotron 
radiation.8,9  The  repetition  rate  of  the  pulses  gives  the  fundamental  frequency, 
whereas  the  pulse  width  determines  the  width  of  the  power  spectrum  band.  The 
power  spectrum  of  mode-locked  lasers  extends  well  above  10  GHz,  an  upper  limit 
in  frequency-domain  spectroscopy  imposed  by  the  optical  detectors  rather  than  the 
light  sources.  The  wavelengths  of  the  above-mentioned  lasers  are  1064  nm  for  the 
NdiYAG,  660-1180  nm  (tunable)  for  the  Tirsapphire,  and  625-780  nm  (tunable) 
for  dye  lasers  using  DCM  or  oxanine  1  dyes.  A  unique  pulsed  source  is  provided 
by  synchrotron  radiation,  which  continuously  covers  the  UV/visible/near-infrared 
spectrum. 
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to-voltage  conversion  and  low-pass  filtering,  the  low-frequency  signal  is  digitized. 
The  A/D  sampling  rate  is  an  integer  multiple  of  Af(nAf),  so  that  the  digital  sam¬ 
pling  is  synchronous  with  the  cross-correlation  signal.  The  sampling  theorem22 
states  that  n  must  be  equal  to  at  least  2  for  the  complete  determination  of  the  har¬ 
monic  at  frequency  A/.  A  higher  sampling  rate  allows  for  the  discrimination  of 
higher  harmonics.  Topical  values  of  n  range  from  4  to  256.  The  digital  signal 
processing  consists  of  a  preliminary  average  of  the  data  sampled  over  a  number 
of  cross-correlation  periods,  followed  by  the  discrimination  of  the  frequency  har¬ 
monic  at  A /  by  discrete  Fourier  transform.  The  first  filter  (averaging)  causes  de¬ 
structive  interference  of  all  the  frequencies  that  are  not  multiples  of  A/,  while  the 
discrete  Fourier  transform  resolves  each  harmonic  of  A/.  The  rejection  of  higher 
harmonics  is  accomplished  by  a  factor  greater  than  2000  (Ref.  [21]).  Figure  7.1 
shows  a  schematic  illustration  of  the  heterodyne  detection  and  digital  acquisition. 

7.2.6  A  FREQUENCY-DOMAIN  TISSUE  SPECTROMETER 

We  now  describe  a  specific  frequency-domain  instrument  for  near-infrared  spec¬ 
troscopy  and  imaging  of  tissues.  After  the  initial  design  and  the  implementation  of 
a  working  prototype  at  the  University  of  Illinois,23  this  instrument  has  undergone 
engineering  and  software  development  at  ISS,  Inc.,  Champaign,  IL.  A  schematic 
diagram  of  the  instrument  is  shown  in  Figure  7.2.  The  light  sources  are  laser  diodes 
that  are  DC  biased  near  threshold  for  laser  operation.  A  superimposed  110-MHz 
radio-frequency  current  signal  modulates  the  optical  output  of  the  laser  diodes 
with  a  modulation  depth  close  to  100%.  The  optical  detectors  are  photomultiplier 


Figure  7.1  Heterodyne  detection  scheme  and  digital  signal  processing  with 
Fourier  filter.  The  intensity  of  the  light  source  is  modulated  at  a  radio-frequency 
(RF)  /,  typically  on  the  order  of  1 00  MHz.  The  beating  between  the  detected  signal 
at  frequency  /  and  the  down-conversion  signal  at  frequency  f  +  Af  determines  a 
low  frequency  component  A /  that  is  filtered,  digitized  (A/D),  and  fast  Fourier  trans¬ 
formed  (FFT)  to  yield  the  amplitude  and  phase  readings.  The  two  radio-frequency 
synthesizers  are  synchronized  (sync). 
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that  the  758-nm  lasers  (now  discontinued)  can  be  replaced  by  690-nm  lasers  (for 
instance,  Hitachi  HL6738MG).  Furthermore,  the  laser  diodes  are  mounted  on  in¬ 
terchangeable  boards,  so  that  one  can  change  the  laser  wavelength  according  to  the 
specific  application.  For  instance,  one  can  use  boards  containing  lasers  emitting 
at  eight  different  wavelengths  for  tissue  spectroscopy  (for  instance,  Figure  7.5  re¬ 
ports  data  collected  at  633,  670, 751,  776,  786,  814,  830,  and  841  nm  (Ref.  [24])). 
Alternatively,  one  can  use  a  laser  board  containing  diodes  emitting  at  the  same 
wavelength  to  maximize  the  number  of  source  locations  in  single-wavelength  opti¬ 
cal  imaging.  In  general,  the  instalment  provides  2  parallel  detectors  and  16  sources, 
resulting  in  32  source-detector  channels.  The  optical  probes  shown  in  Figure  7.2 
are  designed  to  perform  quantitative  tissue  oximetry.  They  implement  the  multi¬ 
distance  measurement  scheme10,25  that  was  found  to  be  robust  and  accurate  for 
in-vivo  applications.  This  scheme  consists  of  placing  multiple  (in  this  case,  four) 
source  fibers  at  different  distances  from  the  single  detector  fiber.  In  this  fashion, 
one  measures  the  spatial  dependence  of  the  DC  intensity,  AC  amplitude,  and  phase 
data,  which  is  related  to  the  tissue  optical  coefficients  by  a  diffusion  model  (see 
Section  7.3.7).  If  two  wavelengths  are  used,  each  illumination  location  consists  of 
two  source  fibers  (one  per  wavelength).  One  of  the  advantages  of  the  multidistance 
scheme  is  its  insensitivity  to  motion  artifacts,  which  allows  for  measurement  while 
the  subject  moves  or  exercises  [see  Figure  7.6c],  Multiple  laser  diodes  time-share 
the  optical  detector  by  means  of  a  multiplexing  circuit  that  turns  the  laser  diodes 
on  and  off  in  sequence,  so  that  only  one  laser  is  on  at  any  given  time.  The  mul¬ 
tiplexing  rate,  which  determines  the  on-time  of  each  laser  diode,  is  adjustable  by 
software.  A  typical  multiplexing  rate  of  100  Hz  corresponds  to  a  10-ms  on-time  per 
diode.  The  laser  rise  time,  determined  by  protection  circuitry,  is  about  1  ms.  There¬ 
fore,  the  data  collected  during  the  first  1  ms  of  the  on-time  period  of  each  laser 
are  discarded.  During  a  typical  10-ms-long  laser  on-time,  45  cross-correlation  pe¬ 
riods  are  processed  (the  5-kHz  cross-correlation  frequency  corresponds  to  a  period 
of  0.2  ms).  As  few  as  25  cross-correlation  periods  (5  ms)  can  be  processed  per 
laser  on-time,  leading  to  a  maximum  multiplexing  rate  of  about  170  Hz.  Faster 
computers  and  higher  cross-correlation  frequencies  may  further  increase  the  data 
acquisition  rate.  A  number  of  multiplexing  cycles  (selected  by  software)  can  be 
averaged  to  increase  the  signal-to-noise  ratio  when  fast  measurements  are  not  re¬ 
quired. 

The  instrumental  noise,  which  depends  on  the  acquisition  time,  is  typically 
much  smaller  than  the  physiological  fluctuations  observed  in  vivo.  The  noise  can 
be  estimated  by  the  standard  deviation  of  the  temporal  fluctuations  of  the  mea¬ 
sured  data.  Figure  7.3  shows  typical  temporal  traces  of  the  average  intensity,  am¬ 
plitude,  phase,  absorption  coefficient,  and  reduced  scattering  coefficient  measured 
with  an  acquisition  time  per  diode  of  160  ms  (average  of  sixteen  10-ms  cycles)  on 
a  tissuelike  solid  phantom  (made  of  silicone).  The  optical  coefficients  of  the  phan¬ 
tom  [Figure  7.3d]  match  typical  optical  properties  of  blood-perfused  tissues  in  the 
near-infrared.  The  standard  deviation  errors  in  the  various  readings  are  reported  in 
Table  7.1.  The  instrumental  errors  decrease  for  longer  acquisition  times. 
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7.3  Modeling  Light  Propagation  in  Scattering  Media 

7.3.1  The  Boltzmann  transport  equation  (BTE) 

The  Boltzmann  transport  equation  (BTE)  is  a  balance  relationship  that  describes 
the  flow  of  particles  in  scattering  and  absorbing  media.  The  propagation  of  light  in 
optically  turbid  media  can  be  modeled  by  the  transport  equation,  where  the  photons 
are  treated  as  the  transported  particles.  If  we  denote  the  angular  photon  density 
with  i/(r,  ft,  f),  which  is  defined  as  the  number  of  photons  per  unit  volume  per  unit 
solid  angle  traveling  in  direction  ft  at  position  r  and  time  t,  we  can  write  the  BTE 
as  follows:26 

— uft  ■  Vw( r,  ft,  t )  —  v([ia  +  Pj)n(r,  ft,  t) 

+  up,  J  u(r,  ft',  f)/(Q,1  ft)dft'  +  q( r,  ft,  t),  (7.1) 

471 

where  v  is  the  speed  of  light  in  the  medium,  \ia  is  the  absorption  coefficient  (units 
of  cm-1),  p5  is  the  scattering  coefficient  (units  of  cm-1),  /(ft',  ft)  is  the  phase 
function  or  the  probability  density  of  scattering  a  photon  that  travels  along  direc¬ 
tion  ft'  into  direction  ft,  and  q(r,  ft,  t)  is  the  source  term.  q( r,  ft,  t)  has  units  of 
s-1  m-3  sr-1  and  represents  the  number  of  photons  injected  by  the  light  source 
per  unit  volume,  per  unit  time,  per  unit  solid  angle  at  position  r,  time  t,  and  di¬ 
rection  ft.  The  left-hand  side  of  Eq.  (7.1)  represents  the  temporal  variation  of  the 
angular  photon  density.  Each  one  of  the  terms  on  the  right-hand  side  represents 
a  specific  contribution  to  this  variation.  The  first  term  is  the  net  gain  of  photons 
at  position  r  and  direction  ft  due  to  the  flow  of  photons.  The  second  term  is  the 
loss  of  photons  at  r  and  ft  as  a  result  of  collisions  (absorption  and  scattering).  The 
third  term  is  the  gain  of  photons  at  r  and  ft  due  to  scattering.  Finally,  the  fourth 
term  is  the  gain  of  photons  due  to  the  light  sources.  Let  us  now  define  some  of  the 
quantities  used  to  describe  photon  transport. 

Angular  photon  density :  w(r,  ft,  t) 

u( r,  ft,  t)  is  defined  such  that  u( r,  ft,  t)drdQ  represents  the  number  of  photons 
in  dr  that  travel  in  a  direction  within  c/ft  around  ft.  The  units  of  u(r,  ft,  t) 
arem-3sr-1. 

Photon  radiance:  L(r,  ft,  t) 

L(r,  ft,  t)  =  vu{ r,  ft,  t).  L(r,  ft,  f)c/ft  represents  the  number  of  photons  traveling 
per  unit  time  per  unit  area  (perpendicular  to  ft)  in  a  range  of  directions  within  c/ft 
around  ft.  The  units  of  L(r,  ft,  f)  are  s-1  m-2  sr-1 . 

Photon  density:  U  (r,  t) 

U (r,  f)  =  f4nu(r,  ft,  t)dQ.  The  photon  density  is  the  number  of  photons  per  unit 
volume.  The  units  are  m-3 . 


du(r,  ft,  t) 
ft 
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nomials,  by  recalling  that  a  function  H(x)  (which  is  sectionally  continuous  to¬ 
gether  with  its  derivative  in  the  interval  [—1, 1])  has  the  general  Legendre  series 
representation31 

00  21  +  1 

H(x)  =  Y/~-HlPl(x),  (7.5) 

/= o 

where  Pi(x )  is  the  Legendre  polynomial  of  order  /,  and  Hi  =  H (x')Pi(x')dx' . 

We  then  write 


00  9  7  4-  1 

/<n'  •  s>  =  ■  m 

h  4,1 

oo  l 

=  E  E  fiYf  &)¥?(&),  (7.6) 

1=0  m=—l 


where  the  last  expression  follows  from  the  addition  theorem  for  spherical 
harmonics,31  namely  Pi(&  ■  £2)  =  4n/(2l  -f  l)YC!m=-i  Y™*  (£2')Y[n(Q).  Here, 
//  =  2n  f*  l  f  (cos  y )  />,  (cos  y  )d  (cos  y ) . 

By  substituting  these  expressions  into  Eq.  (7.1),  we  obtain 


oo  l 

EE 


ir± 


>]* 


+  •  V  +  u(|xfl  +  m)  \uim(r,  t)F,m(£2)  -  qim{r,  t)F/"(£2) 

oo  V 


/=0  m=—l 1  L 

-VM.  f  «/m(r,t)rr(^)E  E  fi'Yf  (&)Y? <&)d& 

■  l'=0  m’=-l' 


=  0.  (7.7) 


The  integral  in  d£l'  can  be  calculated  using  the  orthogonality  property  of  the  spher¬ 
ical  harmonics:  f4n  Yf1*  (£l)d£l  =  8u’8mm>.  The  BTE  thus  becomes 


E  E  (Ia7  +  ,;f2'  v  +  v[^(1  -//)  +  M-fl]]^m(r,O-7/m(r,oV,m(f7)  =  0. 

1=0  m=—l  ' 1  ‘  ' 

(7.8) 

We  then  multiply  this  equation  by  Y^1'  (Cl)  and  integrate  over  Cl  to  obtain  the 
relationship  between  the  specific  coefficients  ulm  and  qLM,  and  all  the  coefficients 
of  the  spherical  harmonic  expansion  of  u ,  mm : 


— mlm(l  0  +  u[|i*(l  -  /,)  +  [ia]uLM( r,  t) 


00  l 


+ 


S2  •  V«/m(r,  t)Y["(n)Y^(^)dtl  =  qLM( r,  t).  (7.9) 


/=0  m=-l 


471 
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— ulm (r,  t)  +  u[Ml  -  fL)  +  V-a]uLM( r,  0 


vuY.hr,  t) 


1  [(L-M+l)(L-M  +  2)-]l/2{  3  .  3  \  M  , 

+  2L  (2L  +  l)(2L+3)  J  \3;c  ldy)VUL+^ 

ir(L  +  M)(£  +  M-l)]1/V  3  .  3  \  M—if  „ 

2[_  (2L  +  1)(2L  -  1)  J  (,3a:  ‘  3}-  J^-1  (r’ 

lfCL  +  M  +  ^CL  +  M+Dl^V  3  ,  .  d  \  M+, 

2[  (2L  +  l)(2L+3)  J  U*  +  *3y/  L+1 

ir(L-M-l)(L-M)l1/2/ 3  .  3  \  M+1/  x 

+  2L  (2L  +  l)(2L  —  l)  J  \3* +*3;y/  1-1  r,? 
UL  +  M  +  l)(L  -  M  +  d_  m 
L  (2L  +  l)(2L+3)  J  dz  1+1 


i<+7(r,  0 


)l1/23 

— i 

3z 


'q-/ 

.(2  L- 


-M)(L  +  M)11/2  3  m  , 


7.3.3  The  p^  approximation 

The  expansion  of  the  BTE  into  spherical  harmonics  has  led  to  an  infinite  set  of 
equations  with  indices  L  (ranging  from  0  to  oo)  and  M  (ranging  from  —  L  to  L). 
Truncation  of  the  Laplace  series  at  L  =  N,  leads  to  the  so-called  approxima¬ 
tion.  The  reason  for  this  name  is  that  the  last  term  in  the  truncated  Laplace  series 
contains  which  can  be  written  in  terms  of  the  associated  Legendre  func¬ 

tions  Pj!f(x),  which  in  turn  can  be  written  in  terms  of  the  Legendre  polynomial 
Pn(x).  The  relationships  are  the  following: 


47t  (N  +  M)\ 


PZ  (cos0)e 


pM(X)=(i-xy/2-fMpNM. 


7.3.4  THE  Pi  APPROXIMATION 


We  now  consider  the  P\  approximation,  which  is  often  used  to  describe  photon 
migration  in  tissues.  In  the  Pi  approximation  ulm( r,  0  is  set  to  0  for  L  >  1.  In 
the  Pi  approximation,  Eq.  (7.13)  is  a  set  of  four  equations.  The  first,  for  L  =  0, 
M  =  0: 


9  1  [2  (3  3  \ 

— M0.o(r,  t)  +  u[p.*(l  -  /o)  +  lia]Mo,o(r,  0  -I-  —  -  i— JuMi  _i(r,  t) 

-\n (s+7)“ii(r-‘,'fyi^”“10(r'')=®0<r,'):  <716> 
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y^(-«u( r,  t)  +  mi  _i(r,  t))\ 

-  #1.1,  t)  +  mi  _i(r,  r))y  +  Mi,0(r,  t) z 


(7.21) 


The  set  of  four  Eqs.  (7. 16)— (7. 19)  of  the  P\  approximation  are  thus  equivalent  to 
the  following  two  equations  (one  scalar  and  one  vectorial): 


■~U (r,  t )  +  u[n*(l  -  /o)  +  \ia]U(r,  0  +  V  •  J(r,  t)  =  VSt^o.o,  t), 
“J(r,  t)  +  [m(l  -  /i)  +  M-a]J(r,  o  +  f) 


2^1-1, 0-?u,  0)x 


-  V  r  (91,-1, 0  +  91,1, 0)y  +  91,0,  t)i 


(7.22) 


(7.23) 


The  vectorial  equation  is  obtained  by  combining  Eqs.  (7.17),  (7.18),  (7.19)  accord¬ 
ing  to  the  following  formal  relationship: 

V^/3[(7.17)  -  (7.19)]x  -  iv^/3[(7.17)  +  (7.19)]y  +  V4^/3(7.18)z. 


From  the  general  definition  of  the  coefficients  /;,  we  find  that  /o  and  f\  are  given 
by 


1  1 

f(cosy)Po(cosy)d(cosy)  =  2nj  f(cosy)d(cosy)  —  l,  (7.24) 


fa  =  27t  j 

-1 

1 

/1  =  2n  j  f(cosY)P\(cosy)d(cosy) 


-1 


-1 


=  Znf 


/(cosy)  cos  yd  (cosy)  =  (cosy), 


(7.25) 


where  in  Eq.  (7.24)  we  have  used  the  fact  that  the  scattering  probability  is  nor¬ 
malized  according  to  the  condition  ■  tl)dQ'  =  1,  which  is  equivalent  to 

27t /_ j  /(cosy)J(cosy)  =  1.  Therefore  /o  is  1,  whereas  f\  is  the  average  cosine 
of  the  scattering  angle  y((cosy)).  The  source  terms  in  Eqs.  (7.22)  and  (7.23)  are 
formally  a  monopole  term  (spherically  symmetric)  and  a  dipole  term,  respectively. 
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The  average  values  of  the  coordinates  of  P3  are32 


(*3)  =  (b)  -  0, 

(23)  =  —  (1  +  <cos 9)  +  (cosd)2). 


(7.31) 


In  general,  at  the  n-th  order  of  scattering,  the  average  values  of  the  coordinates  of 
the  scattering  point  P„  are32 


(xn)  =  (>«)  =0, 


1  "-1 

{Zn)  =  ~ ^(cosy)*  = 
Pj 


k= 0 


1  —  (cosy)" 
Ml  “  (cosy))’ 


(7.32) 


where  we  have  used  the  result  for  the  geometric  series  J2kZ 0  ak  =  (l  —  an)/(l  —  a), 
with  a  <  1.  In  the  limit  of  a  high  number  of  scattering  events  (n  ■— >-  00),  (xqo)  = 
( y 00)  —  0  and  (z0 o)  =  1/[Pj(1  —  (cosy))]  give  the  coordinates  of  the  center 
of  symmetry  of  the  statistical  photon  distribution.  In  particular,  the  coordinate 
(Zoo)  can  be  interpreted  as  the  average  distance  between  consecutive,  effectively 
isotropic  scattering  events,  and  its  inverse  is  defined  as  the  reduced  scattering  co¬ 
efficient  |x' : 


t4  =  H*(l-{cosy>).  (7.33) 

In  the  case  of  isotropic  scattering,  (cos  y )  =  0  and  |x'  =  .  In  the  case  of  forward 

scattering,  (cosy)  =  1  and  p'  =  0. 


7.3.6  THE  Pi  EQUATION  AND  THE  STANDARD  DIFFUSION 
EQUATION  (SDE) 


We  now  reduce  the  Pi  approximation  to  a  single  equation  for  the  photon  density 
U (r,  t).  From  Eq.  (7.27)  we  obtain  J(r,  t): 


J(r,  t )  = 


1 


-J(r,  0 


VU(r,  t)  + 


1 


v{\i's  +  \ka)dt~'"''  3(p'+pfl) 

~  ■DVC/(r>  0  +  — Si(r,  0. 
vl  3 1  v 


(Pi  +  Pa) 


Si(r,  f) 


(7.34) 


where  we  have  defined  the  diffusion  coefficient  D  =  v/[3(p(.  +  pa)].  By  substitut¬ 
ing  this  expression  for  J(r,  t)  in  Eq.  (7.26),  we  get 

d  3  D  9 

-C/(r,  t)  +  v\iaU(r,  t)--y-V.  J(r,  t)  -  D V2C/(r,  t) 
d  t  vz  a t 
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and  the  photon  flux  J(r,  t)  is  related  to  the  photon  density  U (r,  t)  by  Fick’s  law: 

J(r,  t)  =  — DVI7 (r,  t).  (7.40) 

In  the  frequency  domain,  3/3 1  -*■  —ico  and  the  diffusion  equation  takes  the  form 
of  the  Helmholtz  equation: 


(V2  +  k2)U(  r)  =  — (7.41) 

where  k 2  =  (ico  —  v\ia)/D. 

7.3.7  Solution  of  the  standard  diffusion  equation  in  the 

FREQUENCY  DOMAIN 


The  solution  to  the  diffusion  equation  for  a  homogeneous,  infinite  medium  con¬ 
taining  a  harmonically  modulated  point  source  of  power  P  ( co )  at  r  —  0  is  given 
by33 


U  (r,  co)  — 


P(co)  eikr 
AnD  r 


(7.42) 


The  explicit  expressions  for  the  average  photon  density  ((/dc)»  and  for  the  ampli¬ 
tude  (C/ac)  and  phase  (4>)  of  the  photon-density  wave  are25,34,35 


UDC{r)  = 


Pdc 

AnD  r 


{/AC(r,  co)  = 


_r(vtifl/2D)1/2[(]  +  -S^.)l/2+l]l/2 

P(co)e 

AnD  r 


<$>{r,(o)=r{v\ial2D)112  ^1  + 


CO 

v2\i2a 


1/2  -,1/2 

-  1 


+  4*5 


(7.43) 

(7.44) 

(7.45) 


where  <f>s  is  the  source  phase  in  radians.  Analytical  solutions  in  the  frequency  do¬ 
main  have  also  been  reported  for  a  semi-infinite  medium,25’35,36  infinite  slab,37 
cylindrical,  and  spherical  geometries.37  Equations  (7.39)  and  (7.41)  refer  to  homo¬ 
geneous  media.  For  quantitative  tissue  spectroscopy  and  oximetry  (Section  7.4), 
one  typically  assumes  that  tissues  are  macroscopically  homogeneous,  so  that 
Eqs.  (7.39)  and  (7.41)  are  applicable.  By  contrast,  optical  imaging  of  tissues  (Sec¬ 
tion  7.5)  aims  at  measuring  the  spatial  distribution  of  the  tissue  optical  properties, 
and  Eq.  (7.39)  must  be  generalized  to  account  for  the  spatial  dependence  of  \ia 
and  D. 
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Wavelength  (nm) 


Figure  7.4  Absorption  spectra  of  the  three  dominant  near-infrared  chromophores 
in  tissues,  namely  oxy-hemoglobin  (HbC>2),  deoxy-hemoglobin  (Hb),  and  water 
(H2O).  The  absorption  coefficient  is  defined  to  base  e.  The  concentrations  of  Hb 
and  Hb02  are  both  assumed  to  be  50  pM,  a  typical  value  for  blood-perfused  tis¬ 
sues.  These  spectra  are  obtained  from  compiled  absorption  data  for  water42  and 
hemoglobin.43 

by  a  photon  before  being  absorbed.  In  the  NIR,  typical  values  of  pa  in  tissues  range 
from  0.02  to  0.30  cm-1 .  The  photon  mean  free  path  for  absorption  thus  ranges  be¬ 
tween  about  3  and  50  cm. 

7.4. 1.2  Scattering 

The  scattering  properties  are  determined  mainly  by  the  size  of  the  scattering  par¬ 
ticles  relative  to  the  wavelength  of  light,  and  by  the  refractive  index  mismatch 
between  the  scattering  particles  and  the  surrounding  medium.  In  biological  tissues, 
the  scattering  centers  are  cells  and  cellular  organelles.  In  the  medical  spectral  win¬ 
dow  (700-900  nm),  cellular  organelles  have  dimensions  comparable  to  the  wave¬ 
length,  and  their  index  of  refraction  is  relatively  close  to  that  of  the  cytosol  and 
extracellular  fluid.  As  a  result,  light  scattering  in  tissue  is  mainly  forward  directed 
(i.e.,  the  scattering  angle  y  introduced  in  Section  7.3.2  is  less  than  90  deg)  and 
shows  a  weak  wavelength  dependence. 

As  discussed  in  Sections  7.3.1  and  7.3.5,  the  scattering  properties  of  tissues  are 
described  by  two  parameters:  (1)  the  scattering  coefficient  (py),  defined  as  the  in¬ 
verse  of  the  average  photon  path  length  between  successive  scattering  events;  and 
(2)  the  average  cosine  of  the  scattering  angle  (/j  =  (cosy)).  From  the  definition 
of  11^,  it  follows  that  l/ps  is  the  average  distance  traveled  by  a  photon  between 
successive  scattering  events,  as  shown  by  Eq.  (7.28).  Even  though  each  scattering 
event  is  mainly  forward  directed,  after  a  number  of  collisions  a  photon  loses  mem¬ 
ory  of  its  original  direction  of  propagation.  Under  these  conditions,  we  can  say  that 
the  photon  has  experienced  an  effectively  isotropic  scattering  event.  As  discussed 
in  Section  7.3.5,  the  reduced  scattering  coefficient  [p'  =  (1  -  /i)Pj]  represents 


Frequency-Domain  Techniques  for  Tissue  Spectroscopy  and  Imaging 


23 


If  the  extinction  spectra  £;  (A.)  of  all  N  species  are  known,  the  concentrations  Q 
can  be  determined  by  measuring  \ia  at  N  or  more  wavelengths,  so  that  the  lin¬ 
ear  system  of  Eq.  (7.46)  is  fully  determined.  This  approach  requires  that  \xa  be 
measured  independently  of  \i!s .  In  the  frequency  domain,  Eqs.  (7.43)-(7.45),  or  the 
corresponding  equations  for  semi-infinite  or  other  geometries,  can  be  used  to  quan¬ 
titatively  measure  \ia  and  p.'  ,36  The  spatial10  or  the  frequency17  dependence  of  the 
DC  intensity,  AC  amplitude,  and  phase  can  also  be  used  effectively  for  in-  vivo  tis¬ 
sue  spectroscopy.  It  is  worth  noting  that  this  approach  to  tissue  spectroscopy  yields 
quantitative,  absolute  values  of  the  concentration  of  chromophores  in  tissue. 

7.4.3  Near-infrared  tissue  oximetry 

In  many  cases,  the  absorption  spectra  of  tissues  can  be  well  described  by  consid¬ 
ering  only  three  chromophores;  namely,  oxy-hemoglobin,  deoxy-hemoglobin,  and 
water.  For  example,  Figure  7.5  shows  the  absorption  coefficients  of  human  mus¬ 
cle  (forearm)  and  brain  (forehead)  measured  at  eight  NIR  wavelengths  (squares) 
under  resting  conditions.  The  lines  in  Figure  7.5  are  the  best  fit  absorption  spectra 
corresponding  to  a  linear  combination  of  the  water,  oxy-hemoglobin,  and  deoxy- 
hemoglobin  extinction  spectra.  In  the  fits,  the  water  concentration  (by  volume)  is 
assumed  to  be  75%  for  the  muscle,  and  80%  for  the  brain,45  while  the  concentra¬ 
tions  of  oxy-hemoglobin  ([HbC^])  and  deoxy-hemoglobin  ([Hb])  are  the  parame¬ 
ters  for  fit.  The  best  fit  spectra  of  Figure  7.5  yield  values  of  [Hb02]  and  [Hb]  of 


Figure  7.5  Quantitative  absorption  spectroscopy  of  tissues  in  vivo .24'46  The 
squares  represent  the  eight-wavelength  spectra  measured  with  frequency-domain 
spectroscopy  in  the  skeletal  muscle  and  in  the  forehead  (according  to  the  label)  of 
human  subjects.  The  lines  are  best  fits  using  a  linear  combination  of  the  extinction 
spectra  of  oxy-hemoglobin,  deoxy-hemoglobin,  and  water. 
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(a) 


(b) 


Figure  7.6  Temporal  traces  of  optically  measured  total  hemoglobin  concentration 
(THC  =  [HbCh]  +  [Hb])  and  tissue  saturation  (StC>2)  in  human  subjects.  The  mea¬ 
surements  were  conducted  on  a  forearm  muscle  during  (a)  arterial  or  (b)  venous 
occlusion,  and  (c)  on  a  calf  muscle  during  treadmill  exercise.  In  panel  (c),  the  step 
function  represent  the  treadmill  speed,  which  was  changed  by  steps  of  1  mph 
every  2  minutes.50  Running  started  at  a  speed  of  4  mph.  The  frequency-domain 
measurements  allow  for  absolute  y-axis  readings 


Near-infrared  oximetry  lends  itself  to  measurements  on  the  human  head  as 
well.  Figure  7.7a  shows  the  arrangement  of  optical  fibers  on  the  forehead  of  a 
human  subject.  The  cerebral  oxygenation  may  be  varied  by  changing  the  fraction 
of  oxygen  inspired  by  the  subject.  The  resulting  changes  in  the  cerebral  tissue 
oxygenation  (StC>2)  are  reported  in  Figure  7.7b.  Figure  7.7b  also  shows  a  com¬ 
parison  between  two  traces  of  arterial  saturation  (SaC>2);  one  is  measured  with  a 
commercial  pulse  oximeter,  while  the  other  is  measured  with  a  frequency-domain 
tissue  oximeter.46  Both  instruments  use  the  oscillatory  components  of  the  opti¬ 
cal  signals  at  the  frequency  of  the  heartbeat  (which  are  directly  associated  with 
the  arterial  pulsation,  and  therefore  with  arterial  blood  volume).  The  key  differ- 
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Figure  7.8  Frequency-domain,  near-infrared  oximetry  of  the  newborn  infant  brain, 
(a)  Photograph  showing  the  application  of  a  flexible  optical  probe  to  the  infant’s 
head,  (b)  Temporal  traces  of  cerebral  tissue  saturation  (StC>2)  measured  with 
frequency-domain  optical  oximetry,  and  systemic  arterial  saturation  (Sa02)  mea¬ 
sured  with  a  pulse  oximeter  at  an  infant’s  toe.  Deoxygenation  is  achieved  by 
a  controlled  reduction  in  the  fraction  of  inspired  oxygen  (in  collaboration  with 
G.  P.  Donzelli  and  S.  Pratesi,  University  of  Florence,  Italy). 


duced  by  a  controlled  decrease  in  the  fraction  of  inspired  oxygen.  The  arterial 
saturation  and  the  cerebral  saturation  show  a  qualitatively  similar  trend.  However, 
a  local  cerebral  deoxygenation  would  be  detected  only  by  the  NIR  cerebral  mea¬ 
surement.  The  manual  application  of  the  optical  probe  shown  in  Figure  7.8a  can  be 
replaced  by  a  flexible  headband57,58  or  by  an  infant  cap  for  continuous  long-term 
monitoring. 

7.4.4  Measurements  of  optical  scattering  in  tissues 

Traditionally,  optical  studies  have  been  targeted  at  measurements  of  the  absorption 
properties  of  tissues.  The  strong  scattering  of  light  in  tissue  constitutes  a  com¬ 
plication  in  quantitative  absorption  spectroscopy.  However,  it  has  recently  been 
suggested  that  the  reduced  scattering  coefficient  itself  may  provide  information 
about  physiologically  relevant  parameters.  For  instance,  it  has  been  shown  that 
mitochondria  are  the  main  source  of  light  scattering  in  the  liver,  and  possibly  in 
other  tissues  as  well.59  Since  a  number  of  metabolic  processes  related  to  cellular 
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effects.  Theoretical  and  experimental  studies  of  the  optical  properties  of  the  human 
sclera  (the  white  outer  coat  enclosing  the  eyeball)  have  also  shown  a  significant 
scattering  change  induced  by  glucose  and  by  other  osmotically  active  species.63 
While  the  monitoring  of  glucose  concentration  by  NIR  spectroscopy  is,  at  the 
present  stage  of  development,  far  from  being  clinically  applicable,  it  nevertheless 
points  out  the  diversity  of  potential  diagnostic  applications  of  medical  optics. 

7.5  Optical  Imaging  of  Tissues 
7.5.1  General  concepts 

Optical  imaging  relies  on  its  sensitivity  to  optical  properties  of  tissues,  namely  ab¬ 
sorption  and  scattering.  Consequently,  the  contrast  in  NIR  imaging  originates  from 
spatial  variations  in  the  optical  absorption  and  scattering  properties  of  the  tissue. 
These  spatial  variations  can  be  due  to  a  local  change  in  hemoglobin  concentration 
or  saturation  (for  instance,  a  hematoma  or  heterogeneous  vascularization,  blood 
flow,  or  oxygenation),  a  localized  change  in  the  tissue  architecture  (for  example, 
a  result  of  microcalcifications),  or  the  concentration  of  cellular  organelles.  In  any 
case,  NIR  imaging  is  sensitive  to  physical  properties  of  tissues  that  are  different 
than  those  probed  by  current  diagnostic  imaging  modalities  such  as  x-ray  tomog¬ 
raphy  (mass  density,  atomic  number),  ultrasound  (acoustic  impedance),  magnetic 
resonance  imaging  (proton  density,  nuclear  relaxation  times),  and  positron  emis¬ 
sion  tomography  (accumulation  of  a  radioactive  tracer).  In  this  sense,  NIR  imaging 
is  not  necessarily  an  alternative  to  these  imaging  modalities  but  could  constitute  an 
adjunct  technique  that  complements  existing  methods.  We  point  out  here  that  the 
promise  of  optical  tomography  is  not  in  achieving  a  high  spatial  resolution  (which 
is  intrinsically  limited  by  the  diffusive  nature  of  light  propagation  in  tissues),  but 
rather  in  achieving  high  contrast  and  specificity. 

NIR  tissue  imaging  aims  at  generating  spatial  maps  that  display  either  struc¬ 
tural  or  functional  properties  of  tissues.  Since  optical  imaging  studies  the  spatial 
distribution  of  the  optical  properties,  Eq.  (7.39)  must  be  replaced  with  the  diffusion 
equation  for  inhomogeneous  media: 

dU<?'  -  -  V  •  Z9(r)V£7 (r,  t)  +  v\ia (r)C/(r,  t )  =  q( r,  t).  (7.47) 

ot 

Equation  (7.47)  can  be  used  to  solve  the  forward  problem,  which  consists  of  find¬ 
ing  the  photon  density  U  (r,  r)  corresponding  to  a  given  spatial  distribution  of  the 
optical  properties  and  of  the  photon  sources.  Analytical  solutions  of  Eq.  (7.47) 
are  available  only  for  a  few  inhomogeneous  cases  such  as  those  of  spherical33 
and  cylindrical64  inclusions.  For  arbitrary  inhomogeneous  cases,  Eq.  (7.47)  can  be 
solved  using  numerical  methods  such  as  the  finite-difference  method  (FDM)65  or 
the  finite-element  method  (FEM).66  Alternatively,  a  perturbation  expansion  in  \ia 
and  D  leads  to  a  solution  of  Eq.  (7.47)  in  terms  of  a  volume  integral  involving 
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Another  approach,  using  a  hybrid  of  CW  and  frequency-domain  techniques  has 
been  employed  to  perform  real-time  imaging  (image  acquisition  rate:  5.2  Hz)  of  an 
experimental  brain  hemorrhage  in  a  piglet  model.85  Figure  7.10  shows  the  opti¬ 
cal  image  recorded  noninvasively  through  the  intact  scalp  and  skull  immediately 
following  the  subcortical  injection  (depth  of  injection:  1.2  cm)  of  0.5  cc  of  autol¬ 
ogous  blood.  In  Figure  7.10,  a  photograph  of  a  brain  slice  obtained  after  autopsy 
is  superimposed  to  the  near-infrared  image  to  show  the  geometrical  relationship 
between  the  optical  image  and  the  piglet  brain.  The  optical  image  was  obtained  by 
linearly  combining  the  readings  from  10  source-detector  pairs  according  to  a  linear 
back-projection  scheme.85  While  this  scheme  is  computationally  fast,  it  does  not 
lead  to  quantitative  readings,  so  that  the  absorption  scale  in  Figure  7.10  is  quali¬ 
tative  (i.e.,  the  image  gives  reliable  spatially  resolved  information  on  whether  the 
absorption  increases  or  decreases,  but  the  absorption  changes  are  not  accurately 
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Figure  7.10  Optical  absorption  image  of  a  4  x  4  cm2  area  of  the  piglet  brain,  mea¬ 
sured  non-invasively  immediately  after  the  subcortical  injection  of  0.5  cc  of  blood.85 
The  injection  site  is  indicated  by  the  white  circle.  A  picture  of  a  brain  slice  taken 
after  autopsy  is  superimposed  to  the  optical  image  to  show  the  geometrical  rela¬ 
tionship  between  the  imaged  area  and  the  piglet  brain.  One  can  see  the  increased 
optica!  absorption  over  the  area  surrounding  the  injection  site,  which  corresponds 
to  the  subcortical  hemorrhage. 
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<30  s,98  2.5  min,99  several  hours,57,58  and  not-quantified  “slow  data  acquisition 
rate”100  and  “long  measurement  times.”94  Faster  optical  imaging  approaches,  suit¬ 
able  for  noninvasive  functional  brain  imaging,  have  been  reported  recently  with 
image  acquisition  times  of  500  ms101  and  160  ms.79  Noninvasive  optical  studies  of 
the  brain  are  reviewed  in  Ref.  [102], 

Figure  7.11  shows  the  results  of  a  functional  brain  imaging  study  that  used  a 
multichannel  frequency-domain  optical  instrument  (a  modified  version  of  Model 
No.  96208,  ISS,  Inc.,  Champaign,  IL,  shown  in  Figure  7.2).79  The  optically  mea¬ 
sured  maps  of  cerebral  deoxy-hemoglobin  concentration  show  the  activated  cor¬ 
tical  area  during  voluntary  hand  tapping.  Image  reconstruction  was  based  on  a 
linear  superposition  of  the  optical  signals  from  10  source-detector  pairs.  The  fast 
image  acquisition  rate  (6.25  Hz)  leads  to  online,  continuous  optical  mapping  of 
the  cerebral  cortex  during  the  examination.  Figure  7.11  shows  one  optical  im¬ 
age  every  10  measurements  (i.e.,  every  1.6  s)  for  space  reasons.  The  real-time 
movie  of  the  cerebral  activation  recorded  with  near-infrared  imaging  can  be  viewed 
in  Ref.  [79]. 

Currently,  most  optical  functional  studies  of  the  brain  have  been  based  on 
continuous-wave  (cw)  spectroscopy.  Even  Figure  7.11  shows  results  based  on  the 
DC  component  of  the  measured  amplitude-modulated  signal.  The  reason  for  the 
more  common  use  of  cw  data  in  optical  brain  imaging  is  twofold.  (1)  The  study 
of  cerebral  activation  requires  the  measurement  of  changes  in  the  optical  signals 
that  correlate  with  brain  activity.  Therefore,  even  relative  readings,  i.e.,  changes 
from  an  arbitrary  initial  value,  yield  meaningful  information.  (2)  The  absolute 
measurement  capability  afforded  by  time-resolved  methods  requires  the  applica¬ 
tion  of  a  physical  model  for  light  propagation  inside  tissues.  An  accurate  model 
for  the  highly  inhomogeneous  and  geometrically  irregular  case  of  the  human  head 
is  yet  to  be  achieved.  However,  there  are  a  few  exceptions  to  the  use  of  cw  spec¬ 
troscopy  for  optical  studies  of  the  brain.98,99,103'104  These  studies  show  that  even 
without  absolute  measurements,  time-resolved  data  can  be  extremely  valuable  in 
optical  studies  of  the  human  brain.  In  fact,  the  additional  information  content  of 
time-resolved  data  (the  phase  in  the  frequency-domain,  the  time-of-flight  distrib¬ 
ution  in  the  time-domain)  may  lead  to  more  effective  approaches  to  human  brain 
studies,  for  instance  by  affording  the  optimization  of  the  spatial  region  of  sensitiv¬ 
ity,  and  the  maximization  of  the  sensitivity  to  localized  absorption  and  scattering 
changes. 

7.5.3  Optical  mammography 

From  a  practical  point  of  view,  the  human  breast  lends  itself  to  near-infrared  tis¬ 
sue  spectroscopy  and  imaging.  In  fact,  it  is  easily  accessible  and  its  relatively  low 
hemoglobin  concentration  (10-20  uM)105,106  determines  near-infrared  absorption 
coefficients  (~0.02-0.04  cm-1)  that  account  for  a  deep  optical  penetration  depth 
and  the  optical  transmission  through  the  whole  breast.  Furthermore,  there  is  much 
appeal  in  the  possibility  of  detecting  breast  cancer  using  a  noninvasive  method  that 
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does  not  involve  ionizing  radiation.  The  optical  detection  of  breast  cancer  may  rely 
on  the  perturbations  induced  by  cancer  on  the  host  tissue  rather  than  on  a  direct 
detection  of  optical  signatures  associated  with  cancer  cells. 107  These  perturbations 
include  angiogenesis,107,108  alterations  to  the  blood  flow  and  oxygenation,108  and 
fibroblast  proliferation,107  which  may  yield  modifications  to  the  optical  scattering 
and  absorption  properties.  However,  it  is  not  known  whether  these  cancer-induced 
optical  perturbations  can  lead  to  an  effective  optical  approach  to  breast  cancer  de¬ 
tection  over  the  wide  range  of  cases  observed  in  the  clinical  and  screening  practice 
(different  types  and  stages  of  cancer,  various  benign  tumors,  different  levels  of  tis¬ 
sue  heterogeneity,  etc.).  In  vitro  optical  studies  of  normal  and  diseased  breast  tis¬ 
sues  have  not  found  significant  differences  between  their  optical  properties.109'110 
This  result  further  suggests  that  the  optical  contrast  induced  by  cancer  may  result 
only  from  modifications  to  the  vascularization,  oxygenation,  or  tissue  architecture, 
which  are  significantly  altered  for  the  measurements  in  vitro.  In  fact,  in-vivo  mea¬ 
surements  have  indicated  that  cancerous  breast  tissue  has  a  significantly  larger  ab¬ 
sorption  coefficient  and  a  marginally  larger  reduced  scattering  coefficient  with  re¬ 
spect  to  normal  breast  tissue.111-113  Since  it  is  yet  to  be  demonstrated  that  optical 
mammography  can  rely  solely  on  the  intrinsic  optical  contrast  associated  with  can¬ 
cer,  it  has  recently  been  proposed  to  use  dyes  as  optical  contrast  agents.73,114,1 15 

The  first  clinical  tests  of  breast  transillumination  for  diagnostic  purposes  were 
performed  by  Cutler  in  1929.116  This  approach  was  soon  abandoned  because  of 
the  poor  sensitivity  and  specificity  of  the  method.  In  the  1970s  and  1980s,  techni¬ 
cal  developments  led  to  two  new  optical  techniques  called  diaphanography 1 17  and 
lightscanning, 118  which  induced  a  renewed  enthusiasm  for  optical  mammography. 
These  approaches  employed  a  broad  beam  of  visible  and  NIR  continuous-wave 
light  that  illuminated  one  side  of  the  breast.  On  the  opposite  side  of  the  breast,  the 
examiner  visually  inspected  the  light  transmission  pattern  and  used  a  video  cam¬ 
era  for  image  recording.  The  examination  had  to  be  performed  in  a  dark  room, 
and  the  examiner’s  eyes  needed  to  be  dark-adapted.  Despite  some  encouraging 
initial  results,119  several  clinical  studies  conducted  in  the  late  1980s  have  shown 
that  diaphanography  and  lightscanning  are  inferior  to  x-ray  mammography  both 
as  a  screening  and  as  a  clinical  tool.120,121  As  a  result,  medical  acceptance  of  di¬ 
aphanography  and  lightscanning  has  been  subdued. 

The  introduction  of  time-resolved  methods  for  tissue  imaging  has  yielded 
new  instrumental  tools  for  optical  mammography.  Furthermore,  the  application 
of  the  diffusion  equation  [Eq.  (7.47)]  provides  a  mathematical  model  that  al¬ 
lows  for  a  more  rigorous  approach  to  optical  mammography  with  respect  to  the 
mainly  empirical  studies  of  the  past.  As  a  result,  new  time-domain73,112,122-124 
and  frequency-domain18,105,113,125-128  approaches  for  optical  mammography  have 
recently  been  developed.  Siemens  AG,  Medical  Engineering  (Erlangen,  Germany) 
has  recently  designed  and  clinically  tested  a  frequency-domain  prototype  for  op¬ 
tical  mammography.129  Figure  7.12  shows  a  block  diagram  of  this  prototype  and 
two  representative  optical  mammograms  obtained  on  a  patient  affected  by  breast 
cancer.  The  optical  mammograms  shown  in  Figure  7.12  are  2D  projection  images 
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by  minimizing  the  effects  of  the  breast  geometry  on  the  optical  data.80  As  in  x-ray 
mammography,  the  breast  is  imaged  in  two  directions,  leading  to  a  craniocaudal 
(cc)  and  an  oblique  (ob)  view.  Figure  7.12  shows  that  breast  cancer  can  be  de¬ 
tected  with  excellent  contrast  using  optical  methods.  However,  optical  mammog¬ 
raphy  still  needs  to  improve  its  performance  in  the  detection  of  small  tumors  and 
in  the  discrimination  between  malignant  tumors  and  benign  breast  lesions.  A  pre¬ 
liminary  analysis  based  on  the  criterion  that  an  optica]  mammogram  is  positive  if  it 
shows  a  region  of  abnormal  absorbance  in  both  the  craniocaudal  and  oblique  views, 
has  led  to  a  sensitivity  (fraction  of  cancerous  breasts  successfully  detected)  of  72% 
and  a  specificity  (fraction  of  noncancerous  breasts  correctly  evaluated  as  nega¬ 
tive)  of  52%  on  a  clinical  population  of  131  patients.130  This  result  is  consistent 
with  the  sensitivity  of  73%  obtained  on  69  patients131  with  a  similar  frequency- 
domain  prototype  developed  independently  by  the  research  laboratories  of  Carl 
Zeiss  (Oberkochen,  Germany).132  It  must  be  stressed  that  these  results,  which  are 
based  on  optical  mammograms  like  the  ones  shown  in  Figure  7.12b,  rely  on  2D 
projection  data  at  a  single  wavelength.  The  potential  of  optical  mammography  goes 
beyond  these  capabilities,  so  that  there  are  a  number  of  possibilities  for  developing 
more  effective  approaches.  For  instance,  testing  nonplanar  geometries  without  (or 
with  minimal)  breast  compression,  exploiting  the  spectral  information,  quantify¬ 
ing  and  discriminating  the  breast  absorption  and  scattering  coefficients,  achieving 
a  true  2D  or  3D  spatial  reconstruction  of  the  breast  optical  properties,  performing 
oxygenation  or  blood  flow  imaging,  and  developing  imaging  approaches  based  on 
dynamic  perturbations  (such  as  periodic  light  compression  and  decompression  of 
the  breast),  are  some  of  the  areas  currently  under  investigation  by  research  groups 
worldwide. 

7.6  Future  Directions 

The  key  difference  between  continuous-wave  and  frequency-domain  spectroscopy 
is  the  phase  measurement  capability  afforded  by  the  frequency-domain.  For  optical 
studies  of  tissues,  the  phase  information  can  be  exploited  in  at  least  two  ways. 

First,  by  combining  the  phase  measurement  with  the  DC  intensity,  AC  ampli¬ 
tude,  or  modulation  measurement,  one  can  separately  measure  the  absorption  and 
the  reduced  scattering  coefficients  of  tissues.  This  is  of  paramount  importance  for 
quantitative  tissue  spectroscopy,  imaging,  and  oximetry.  In.  this  chapter,  we  have 
illustrated  the  capability  of  frequency-domain  spectroscopy  to  quantify  the  absorp¬ 
tion  and  the  reduced  scattering  coefficients  [Figures  7.3d,  7.5,  and  7.9],  and  to 
perform  absolute  tissue  oximetry  (Figures  7.6-7. 8). 

The  second  relevant  aspect  of  the  phase  measurement  is  that  the  phase  is  sensi¬ 
tive  to  a  different  tissue  volume  with  respect  to  the  DC  intensity  and  AC  amplitude. 
The  region  of  sensitivity  can  be  evaluated  by  the  effect  on  the  measured  quanti¬ 
ties  of  a  small  perfectly  absorbing  object  located  at  r.  In  particular,  the  DC,  AC, 
and  phase  sensitivities  can  be  expressed  by  the  parameters  eoc(r)  =  DC(r)/DCo, 
eAc(r)  =  AC(r)/ACo,  and  ep(r)  =  O(r)  -  <t>o,  respectively,  where  DC0,  ACq,  Oq, 
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the  figure  caption).133  Figure  7.13  shows  that  the  DC  and  AC  measurement  are 
most  sensitive  to  the  region  along  the  line  joining  the  source  and  the  detector, 
which,  by  contrast,  is  not  the  highest  region  of  sensitivity  for  the  phase.  In  a  non- 
invasive  tissue  measurement,  which  is  better  modeled  by  a  semi-infinite  geometry, 
where  the  illuminating  and  collecting  optical  fibers  are  located  on  the  tissue/air 
interface,  the  phase  measurement  is  typically  less  sensitive  to  the  superficial  tis¬ 
sue  layer  with  respect  to  the  DC  and  AC  measurements.  As  a  result,  the  phase 
information  may  contribute  to  achieving  depth  discrimination  in  optical  imaging 
of  tissues. 

Future  applications  of  frequency-domain  techniques  to  the  optical  study  of  bi¬ 
ological  tissues  will  exploit  the  two  above-mentioned  features,  with  the  objective 
of  developing  innovative  approaches  for  medical  imaging  and  diagnosis. 
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